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ABSTRACT: The neutral B mesons, B and B^ , can oscillate between their particle and 
antiparticle states owing to flavor-changing weak interactions. In recent years, techniques to 
detect these oscillations as a function of the meson's decay time have been developed. In this 
article the physics of flavor oscillations is reviewed and theoretical predictions are summarized. 
The many observations that demonstrate the time-dependence of B'^ — B'^ oscillations are pre- 
sented along with a combined measurement of its frequency, Anid = 0.484 ± 0.015 ps~^. The 
attempts to measure the B^ oscillation frequency, both directly and indirectly, are then summa- 
rized, currently resulting in a limit of Arris > 14.6 ps~^ (95% CL). Finally, values for the CKM 
elements \Vtd\ = (3.6 ± 0.4) x 10"^ and \Vu/Vtd\ > 4.7 (95% CL) are extracted. 
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1 INTRODUCTION 

The ability of a very few neutral mesons to change from their particle to their 
antiparticle state is a remarkable consequence of basic quantum mechanics and 
the structure of the weak interaction. This oscillation from matter to antimatter 
can be used to measure fundamental parameters of the standard model; in addi- 
tion, it might have far reaching effects, such as breaking the matter/antimatter 
symmetry of the universe. 

In 1955, Gell-Mann &: Pais ^ showed that if a meson with strangeness 
5'=— 1, as well as the meson with strangeness S=l exists, then a quantum- 
mechanical mixing due to — K interactions takes place. It produces two 
physical particles, and K2-: that are a mixture of these states of well-defined 
strangeness, or flavor. Gell-Mann & Pais predicted the existence of a long-lived 
neutral kaon, /Cg, which decays to three pions, as the companion to the shorter- 
lived that had already been observed. Lande [||] at Brookhaven confirmed 
the existence of the longer lived state, K^, in 1956. 

If the physical particles are a mixture of states of well-defined flavor, then these 
flavor eigenstates can be considered mixtures of the physical particles. These 
physical states must have slightly different masses, as discussed in the following 
section, and so they develop a phase difference as they evolve in time. Therefore 
the physical particle content of a flavor eigenstate evolves with time - an initially 
pure flavor eigenstate develops a component of the opposite flavor. The mixing 
of flavor eigenstates to form the physical particles, then, is equivalent to the 
oscillations of flavor eigenstates into one another. 

The only hadrons that can undergo these oscillations are the following mesons: 
/C", , and B^. The tt^ is its own antiparticle, the top quark is so heavy that 
it decays before forming stable hadrons, and excited meson states decay strongly 
or electromagnetically before any mixing can occur. 

Such particle- antiparticle mixing has since been seen for B mesons, first in an 
admixture of B^ and B^ by UAl and then in mesons by ARGUS @ and 
later CLEO B^ mixing was established by comparing the time-integrated 
oscillation probability for B^ mesons, measured by ARGUS and CLEO, to that 
measured at LEP, which contains both B^ and B^ contributions Q. Mixing is 
expected to be a very small effect in mesons and has not been observed. The 
experimental part of this review focuses specifically on recent analyses that at- 
tempt to measure the time-structure of B^ and B^ oscillations directly, rather 
than those that studied the oscillations in a time-integrated manner. The latter 
include those from ARGUS and CLEO [Q] and earlier analyses from the experi- 
ments presented in this paper. See for a list of these measurements. 

The outline of this paper is as follows. Section^ collects the various elements 
required to derive the mixing formalism; Section presents results on quantities 
needed to extract Cabbibo-Kobayashi-Maskawa (CKM) matrix elements from 
mixing measurements; Section |3| presents an overview of experimental techniques 
used to measure the time dependence of the flavor oscillations; Section ^ presents 
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measurements of the oscillation frequency; Section ^ presents limits on direct 
searches for B^ flavor oscillations; and Section ^ deals with measuring the Bg 
oscillation frequency through lifetime differences in the B^ system. 

Since this paper's focus is i?-meson mixing, the formalism is presented in terms 
of b quarks coupled with d, s quarks, though it was first derived for the kaon 
system. B is used generically to mean either B^ or B^ in cases where the result 
applies to both; however for clarity and brevity only the B^ Feynman diagrams 
are shown. Also, the charge conjugate of listed decay modes is implied unless 
explicitly stated otherwise. For convenience we have set h = c = 1 and supress 
the c in the units of momentum and mass. 



2 MIXING FORMALISM 



Particle-antiparticle oscillations are possible because of the flavor-changing term 
of the standard model Lagrangian, 

where Vckm is the CKM matrix |^]. A popular parameterization of this matrix is 
that of Wolfenstein Q, which expands each term in powers of the Cabibbo angle 
A 0.22, shown below to ©(A^): 

Vud Vus Vub\ ( A A\\p-ir,)\ 

Vcd Vcs Vcb \ ^ \ -A(l + iA^XS) 1 - AX^ 
Vtd Vts Vtb J \ AX^il- p-i7]) -AX^ 1 / 

This term in the Lagrangian engenders box diagrams involving an internal loop 
with two W bosons, of the sort shown in Figure |2|, which result in nonzero 
transition matrix elements between B^ — B^ and Bg — B^. 

Before calculating these transition amplitudes within the standard model, the 
relevant general formalism is developed, i.e. the quantum mechanics of a two- 
state, particle-antiparticle system weakly coupled to a continuum. 



2.1 Introduction to Mixing Formalism 

The oscillation effect follows from a simple perturbative solution to Schrodinger's 
equation. Let Tio be the Hamiltonian of the strong interaction. If this were 
the only force, there would be stable states, \B >= \bq > and \B >= \bq > 
{q = d,s), that are eigenvectors of TCq. By the CPT theorem, the masses are 
equal, ms = m-^ = mo, and the Hamiltonian is 

Wo = f 7 ' 

V mo 

When the weak interaction TLw is added, the simple two-state system becomes 
much more complicated, as shown in Figure |l[ The weak force is responsible for 
nonzero matrix elements between the two states to a continuum of states (i.e. 
possible decay modes), which can be split into three groups: (a) states accessible 
only to \B>, i.e. states \a > with <a\Ti.w\B >= 0; (b) states accessible only to 
\B>; and (c) states coupled to both \B> and \B>. In addition, the two discrete 
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states are connected by a direct matrix element W12 =< B\Hw\l^ > or via off- 
shell (i.e. with E 7^ mo) continuum states accessible to both. The Hamiltonian 
H = Ho + Hw has infinite dimensions. 



<BIWIB> 




Figure 1: Schematic of a two-state system coupled weakly to a continuum (indi- 
cated by the wiggles) 



The time evolution of a wavefunction that is a pure |i? > at i = is eas- 
ily calculated using standard first-order perturbation theory, and in general the 
wavefunction develops both \B > and continuum components. Typically, one 
limits the state space to that spanned by the \B > and \B >, however, and de- 
rives a matrix Ti. (still referred to as a Hamiltonian even though it is clearly not 
hermitian) which allows one to use the Schrodinger equation in this subspace. It 
is easy to show it has the form: 

^ / mo + dE Wi2 + dE^2 \ _i{ r ri2 \ . . 

V + SEt^ mo + SE J 2 V ^ J' ^ ' 

The introduction of the continuum states has several effects: 

1. It introduces a slight mass shift 6E in both \B > and \B >, which can be 
absorbed into the mass M = mo + SE. 

2. It introduces off-diagonal elements M12 = W12 + SE12 in the real part of H 

due to the coupling of the particle and antiparticle via TCw and via off-shell 
continuum states accessible to both \B> and \B> {5E12). 

3. It introduces an imaginary part to H, with diagonal elements given by the 
matrix elements oi \B > and |5 > to the on-shell continuum states unique 
to each. 

4. It introduces off-diagonal elements into the imaginary part of TC given 
by the matrix elements of |i? > and |i? > to the on-shell continuum states 
common to both. 
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Diagonalizing this 7i gives us the eigenvectors and eigenvalues. The eigenstates 
are conventionally denoted H, L, for heavy and light, rather than short and long 
as customary for the kaon system. If 



\Bl> 
\Bh> 



p\B> +q\B> 
p\B> -q\B>, 



and 



then 



P 



2^ 12 



and the eigenvalues are 

^H,L = 'mH,L 
where the masses and widths of these states are 



tI 12 



(2) 



(3) 



(4) 



mH,mL 



which satisfy 



M±Re\ \Mi2\ 



12 



iReiMuTl^) = M± Am/2, 



r±2/mW|Mi2P 



121 



iRe{Mi2Tl^) = r± Ar/2, (5) 



Am'' 



. AT^ 
4 

AmAr 



4(|M: 



12 



121 



and 



4Re{Mi2T*u) 



(6) 



Equation § can be rearranged to give the flavor eigenstate |i? > in terms 
of \B}{ >,\Bl >■ The standard derivation uses their simple time dependence, 
[\BH{t) >= e"*'^^*!^// >], to obtain the time evolution of a state which was a 
pure |i?> at t = 0, expressed in terms of the flavor eigenstates. (This procedure 
implies the existence of a reference frame with both the \Bl> and \Bh > com- 
ponents of the \B > at rest, and such a frame does not, strictly speaking, exist. 
However, setting the relative velocity of the \Bl > and \B}{ > states to zero is an 
excellent approximation. A comment on its validity and the appropriateness of 
describing the system's evolution as a function of time can be found below.) 



\B{t)> 



l-i\BUt)> +\BH{t)>) 



-iMt 



e 2 



(( 



AT, ,Am, 

e ^ e + e * e 



H — (e 4 "e" 2 
P 



6 4 6 



\B> 



\B> 



(7) 



Similarly, the time evolution of a state which is pure |i3 > as i = is 
1 

Yq 

A 



\B{t)> 



[\BL{t)> -\BH{t)> 

iMt 



"e 2' 



J. ■ iA//t i Ar J. -Ami, , 

e 4 *e*~* - e-—^e-'—^) B > 



, AE+ .■ Am i AT , .■ Am J. , 

+ (e— *e*— * + e"— *e-*— > 



(8) 
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The normalizations of these states are given by 



< B{t)\B{t) > dt = - 



1 + \q/p\^ ^ 1 



r2_ Ar2/4 



r2 + Am2 



< B{t)\B{t) > dt 



r2-Ar2/4 r2 + Am2 



(9) 



Let Pm{t) denote the probability that a particle produced as a S oscillated 
(mixed) and decayed B. [i.e. P^it) = ^1 <B\B{t)> |2.] Let P^(t) denote 
the conjugate probability that this particle did not oscillate, that is, it remained 
unmixed (with similar definitions for initial B states) . Then Equations |^ |^ and 
give the following: 



Puit) 



-rt 



Ar^/4 + r^+A 



Ar 

(cosh —^t + cos Amt), 



-rt 



i+k/pP , i-k/pP 
r2-ArV4 r2+A 



(cosh -^t 



cos Amt), 



(10) 



^i'W = _/ .l^lf'' (cosh ^t + cosAmt), 



^r^-ArV4 r^+A 



^™(*) = i+|,/pP i r^7^ (cosh^t-cosAmt). (11) 
-■^ Vr2~ArV4 r2+Am2 J 

Note that these expressions are not symmetric between B and B states. 

These formulae have two limiting cases: neglecting CP violation in the mixing, 
and neglecting the lifetime difference AT (which also in general implies there is 
no CP violation in the mixing). 

Equation can be rewritten as (similarly for the \Bh >) 



\Bl> 



p + q 



{\B>+\B>) + ]^{\B>-\B>] 



and thus {l — q/p)/{l + q/p) = es is a measure of the amount by which | Bh > and 
\Bl> differ from CP eigenstates. is expected to be very small in the standard 
model, 0(10"'^). The current world average for the B^ system is Re(eB)=0.002ib 
0.007 No measurement or limit exists for the B^ system. The limit of no 
CP violation in mixing is thus q/p = 1. In this limit the B and B symmetry is 
regained, and we obtain unmixed and mixed decay probabilities for both B and 
;B of: 

1 / Ar2\ Ar 

Pu,m{t) = -re-^Ml-^j (cosh— 1± COS Ami), (12) 

where the + sign corresponds to Pu- This form is appropriate for B^ mesons, 
which are not expected to be subject to large CP-violating effects. 

On the other hand, even in the presence of CP violation, a simple form can be 
obtained. The lifetime difference between the heavy and light states is expected 
to be small, AT/T < 1% for the 5° and perhaps as large as 25% for the 5° fTO]. 
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From Equation 0, AF = in general only if F12 = 0. In this case, 



p \Mu- iFi2 V M12 



thus \q/p\ = 1. In this AF = limit, the time evolutions from Equations ^ and 
^ become 

\Bit)> = e-'^'e-i' (^cos^t\B > +ie-"f>sm^t\B > 
\B{t)> = e-'^''e~i' (^cos ^t\B>+ie+''l' sin ^t\B>y (13) 

The mixed and unmixed decay probabilities again become equal for the B and 
B mesons: 

Pu,m{t) = ^re-r*(l±cosAmt). (14) 

This form is expected to be appropriate for B^ mesons, for which a large phase 
<j) (the source of mixing-induced CP violation) is possible. 



Equations 12 and 14 will be the expressions used throughout this paper. Their 
time-integrated versions, expressing the probability that a B decays as a B, are 
(with X = Am/F) 

Jo 2 l + x^ 

and in the AF = limit, 

1 

X = -— (16) 

As has been pointed out in the context of neutrino oscillation experiments [0] , 
the oscillations observed by any experiment are oscillations in space not in time. 
That is, one has a source creating a pure B 01 B meson, which may have oscillated 
by the time it reaches a distant detector. In this spatial picture, we have a source, 
very small compared to the oscillation wavelength, which emits a pure B meson. 
The boundary condition that must be imposed, then, is that the probability of 
finding a B meson at the source must vanish for all time, otherwise a pure B 
would not be emanating. The \Bi > and \Bh > components propagate with 
phase e*(^^'^*~P^'^^\ where x denotes the direction of motion. At the origin, the 
only way to ensure the wavefuntion does not change the relative \Bl > —\Bh > 
phase and develop a B component is the condition El = Eh- That is, the B 
meson has a definite energy. The components \Bl> and \ Bh > will have the same 

energy, but different momenta pl,h = ^J^"^ ~ "^i h respectively. This induces 

spatial oscillations that go as e*^^^"^^)^. Thus, if the light and heavy components 
have the same energy but different momenta, they must be moving with different 
velocities, and there is no frame in which both the \Bl > and \Bh > are at rest. 
In the derivation above, however, we put El^h = ^l,h for both states. 

At the time t, the centers of the \Bl> and \Bh > wave packets have separated 
by a distance d = {vl — VH)t = {pi — PH)t/E. As Am is very small compared to 
the mean mass M (3 x lO"^^ GeV vs. 5.3 GeV), this distance is d w "^^jf^ ct, 
where = E"^ — M^. Setting t as large as 10 ps (the B lifetime is approximately 
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1.5 ps), and noting that the typical i?-meson energy in the current experiments 
is 10 — 30 GeV, we see that the wavepacket separation is only 0{W~^^ fim). This 
justifies treating the system as moving with a single group velocity given by the 
mean velocity v = P/E. Thus the spatial dependence of the oscillations, {pi — 
Ph)x, can be converted to a time dependence using j; = vt = P/E [pi+pu) /2E 
to obtain {p\ — p'\j)t/2E = {m\ — m'\j)t/2E. If we substitute E = 7M, and the 
proper time r = t/7, the time dependence goes as (m| — m|^)r/2M = Amr, 
recovering the form used in the standard treatment above. 



2.2 Mixing in the Standard Model 

The time-evolution and decay probabilities presented in the previous section are 
valid for any perturbation to the strong Hamiltonian. This section evaluates the 
matrix elements M12 and ri2, in the context of the standard model. 



h- 
W 
d- 



Vqb {u,C,t) V*ji 



W )W 

{u,c,t) 



v; 



Vn 



b- 

d- 



w 



Qd "Qb 

(a) 



VQd 



[u,C,t) 

b 



qh 



Figure 2: Lowest order box diagrams responsible for — oscillations. 



To lowest order, the matrix element coupling a B^ and B^ is given by the 
Feynman diagrams shown in Figure ^a,6, with similar diagrams for the B^. With 
three generations, the internal quark lines Q, q can be u, c, or t. Of course we must 
take into account that the quark lines coming in and going out are not free - the 
h and light quark d or s are bound into a hadron. These diagrams are evaluated 
using a method first employed by Gaillard & Lee ||l^ in 1974 for the system. 
The Hamiltonian corresponding to the box diagrams is sandwiched between the 
bound-state hadrons, Ai =< B\Hw\B >, (so that M12 = M./{2mB), including 
the normalization factor) , the momenta and masses of the incoming and outgoing 
quarks are set to zero (since they are small compared with Mw), and the internal 
loop momentum integral is performed. 

Some care must be taken in evaluating these diagrams. A good choice of gauge 
to perform the calculation is the generalized renormalizable gauge, R^, with ^ 
finite, introduced by Fujikawa et al |13|. In the standard model, there is an 
unphysical charged scalar (p^, which is a remnant of the charged Higgs that is 
absorbed to give the W boson its mass. In fact, such particles will generally exist 
in any model which generates masses for the vector bosons through spontaneous 
symmetry breaking. These scalars couple to quarks with a strength proportional 
to the masses. For an incoming d-type quark of generation j and outgoing it-type 
of generation i, the vertex factor is ^^^^^ [^dj(l + 75) - rUu^il - IbWu.dj- 
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Vqb q Vqd Vqb q Vqd 



h- 



)w w )w 

Q I 



(a) (b) 
Vqb q V*^ Vqb q Vqd 



\<t> iw \<i> \<t> 

I Q I 



^Qd ^Qb Vqd ^Qb 

(c) (d) 

Figure 3: Examples of box diagrams for 6 — 6 oscillations, including unphysical 
scalar contributions. 



The W and </> propagators in the gauge are 



M^ + ie 



fc2 



fc2 - ^M^ + ie 



While the unitary gauge, Roo, is convenient for evaluating tree-level diagrams 
(the unphysical scalars drop out), in this gauge there are extraneous singular- 
ities generated in the box diagram from the W^-propagator terms, which go as 
kfj^ky/My^. In fact, the Green's functions are unrenormalizable, though one can 
show that the full S-matrix is finite [13|, so extreme care must be taken in per- 
forming the ^ — > oo limit. 

In the discussion below we work in the 't Hooft-Feynman gauge, ^ = 1, in 
which the propagators are particularly simple. Note that the relative minus sign 
in the propagators causes the diagrams in Figures c to add to, rather than 
cancel, the diagrams in Figures ^a, d. 

For example, the basic matrix element for — B^ transitions through the 
diagram in Figure Sa, involving internal quarks of type Q and q is 



iM 



<B 



Oi9 



'64 
1 



(fk 

(2vr)4 



d(p3)7^(l 
1 



1 



75 j 



75 J 



m, 



1 



{k-pi) 



M^ik 



P3> 



7a*(1 -75)ft(P4 



7^(l-75)6(pi) 



Here (pi, P2) are the incoming (6, d) momenta; (ps, ^4) are the outgoing (d, 6) 
momenta, and = Vq^Vq^ and = V*^Vqh are the CKM factors. 

We now invoke the approximation that the external 6 and d quark momenta 
and masses are zero, pi = 0. After a reduction of the Dirac matrices, the matrix 
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iM 



WW 



fc2 



(27r)4 (A;2 - m|)(A:2 - m2)(/fc2 _ m2,)2 ' 



Applying the same assumptions to the diagrams with the scalars, one finds that 
all diagrams reduce to having a common operator structure, Oil = (c?l7'^^l)^) 
multiplied by a convergent integral over the loop momentum. The result is 



MAB{Xq,XQ) = 

with Xq = m2/m^ and 

EwwiXq,XQ) = 
Ew<j,{Xq,XQ) = 

E^cl>{Xq,XQ) = 



9 



1 1 



2 Xq 
1 XqXQ 



EAB{xq,XQ) < B^'lidL^niflBO > 



xl log Xq x| log XQ ^ 1 



1 



2 Xn 



(l-Xq)^ (1-Xq)2 1 
Xq log Xq _ Xq log Xq 

(1 - xqr 



Xn 



+ 



^-XQ 
1 



(1 - Xq)2 1- Xq I - XQ 



E^w{Xq,XQ) 



XqXQ 



EwwiXq,XQ). 



Table |l| provides values of these functions for 
mi = 170 GeV. 



0, 



1.5 GeV, and 



Table 1: Contributions of the W-W, T^-scalar and scalar-scalar diagrams with 
all combinations of internal quark lines. 



Diagram 


u, u 


u, c u, t 


c, c 


c, t 


t, t 


(a) W-W 


1 


1.9952 0.5353 


0.9955 


0.5352 


0.1339 


(b)+(c) W-(l) 








1.45 X 10"^ 


1.02 X 10"^ 


1.196 


(d) 0-,/. 








3.02 X 10"^ 


2.08 X 10-^ 


0.6690 



Because the scalar couples to the quarks with a strength proportional to 
niq/Mw, only the diagrams with two top quarks contribute substantially to the 
diagrams in Figures ^b, c, d. Although the W — W diagrams with u quarks and c 
quarks have large matrix elements, they very nearly cancel when combined with 
their CKM factors because both the u and c quark are light compared to the W . 
Traditionally, the unitarity of the CKM matrix is assumed, and is used to express 
the n-quark terms VudVub — ~^cdV*b ~ ^tdVtb- example, it is not the case 
that the ut matrix element is small, but rather its contribution almost exactly 
cancels the ct contribution via this GIM-like mechanism |14| if the CKM matrix 
is unitary. 

Thus the full matrix element, summing over all internal quark types, becomes 
simply a sum over c and t quarks and different bosons: 



q,Q=c,t 
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^ [EAB{Xq,XQ) 

A,B=W,(h 



EAB{Xg,0) - EABiO,XQ)+EABiO,0)] . 



To summarize, we have made the following assumptions: (a) The W mass is 
generated by a spontaneous symmetry-breaking mechanism, resulting in scalar 
remnants; (6) the momenta and masses of the b and d quarks inside a bound 
meson can be neglected; and (c) the CKM matrix is unitary. 

This program was first carried out, with the important scalar terms included, 
by Inami & Lim |15| in 1981. The result of the sum over bosons is 



f3iXQ,Xg) 



XQXq 



In^-^ 



, InxQ 3 xllnxq 



XQ 

3 



4(1-xq)2 4(1 



Xn 



XQXq 



A{l-XQ){1-Xqy 

For the same quark type on the two internal lines (g 



f2{Xq) 




1 



Q), this reduces to 

3 log Xq 



2(1 



2 1 



\3- 



(17) 



(18) 



Table § shows the size of the Inami-Lim functions f2{xc), f2{xt) and f'i{xt,Xc) 
(= fsixc, Xt)), and the CKM factors that multiply each term. For the cases of 
S° and 5° mesons, it turns out that these CKM factors are almost constant, 
independent of the type of quarks in the internal loop. For reference, the CKM 
coefficients for the system are shown, which vary by five orders of magnitude. 



Table 2: Factors entering the matrix element, which is proportional to the product 
of the Inami-Lim function and the CKM term. 



Internal 
quarks 


I-L factor 


CKM 


S° CKM 


K° CKM 


c,c 


3.5 X 10"'' 




A^A^ 


A^ 






(7.4 X 10-5) 


(1.4 X 10-3) 


(2.7 X 10-2) 


c,t 


3.0 X 10-3 


A'^X^ll- p-irj\ 


A^A^ 


A'^X^ll- p-irj\ 






(7.3 X 10-5) 


(1.5 X 10-3) 


(8.8 X 10-6) 


t,t 


2.5 


A'^X^ll- p-i7]\^ 


A^X^ 


A^X^°\1- p-i7]\'^ 






(7.2 X 10-5) 


(1.5 X 10-3) 


(1.1 X lO-'^) 



Note that the function f2ixq), for which both internal quarks are of the same 
type, is proportional to Xg, i.e. the square of the mass of the internal quark, 
and so the top quark loop completely dominates the charm loop's contribution. 
The function /^{xtjXc) is proportional to Xc and is relatively small, an ~ 0.1% 
correction. 

In the approximation that only the top quark contributes, the matrix element 
can thus be written (with a factor for all permutations) as 

^ = Po\VtdVtb\'mU2{m^/Ml)<B\dL7''bL)W>- (19) 
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Hence, we need to evaluate the effective matrix element 

<B\(dr (I-75) bf\B>= 

J2 < 5^7^(1 - 75)b\n X n|d7^(l - j5)b\B > 

n 

= Bb< B\dY{l - 75)5|0 X 0|d7^(l - 75)6|S > 

= Bb\ < 0|d7''75&|5 > |^ (20) 

where Bb is the bag factor, which is 1 if the vacuum insertion |0 >< 0| saturates 
the sum over all intermediate states. 

Substituting < 0\d^'^^5b\B >= i/bp'^, where Jb is the B decay constant into 
this expression, and taking into account the three colours of the quarks gives 

<B\{dj'^{l-j,)bf\B> = ^Bb/Wb- (21) 

Finally, the diagrams in Figure ^ can have arbitrary numbers of gluons running 
between any of the quark lines. The effect of this on M is encoded in the QCD 
correction factor r]B, which multiplies the A4 above. 

The size of this correction, including next-to-leading-order (NLO) diagrams, 



has been estimated by Buras et al |16| 



rjB = 0.55 ±0.01. (22) 

These authors found that rjB depends upon the definition of the top quark mass 
and that only the product r]Bf2{xt) is insensitive to this choice. To be consistent 
with this estimate of ry^, the top quark mass must be chosen to be mt(mj'°^'^), the 
mass calculated using the MS renormalization scheme evaluated at the measured 
(pole) mass of 174ib5 GeV This gives a mass 6-7 GeV lower than the measured 
mass. Some confusion over rjB exists in the literature; however ry^, by definition, 
is the same for B^ and B^ mesons. 



Putting together Equations 19, 21, 22 results in 



|Mi2| = ^ = ^mBflBBVEmlUml/M^WCM''- 

To relate this to the heavy- light mass difference Am, we in principle also need 
the absorptive part of the matrix element, ri2, as per Equation ^. However, 
ri2 ^ M12 as argued below, and can be neglected. 

The matrix element ri2 involves decays to modes common to both \B > and 
\B > that are on-shell and thus energetically allowed, meaning the top quark 
loops do not contribute. Since the matrix element is given approximately by 
the square of the mass of the state in the loop, and this is « for on-shell 
transitions, we obtain ri2/Mi2 ~ rn^/m^ <C 1. 

Alternatively, this can be seen dirctly from the data. The fractional lifetime 
difference AF/F is expected to be small, as previously mentioned - less than 0.01 
for the B^ and up to 0.25 for the Bg. As shown in Sections |^ and ^, current 
experiments demonstrate that Am ~ 0.7F for the B^ system, and Am^ > 20F 
for the Bg. Combining this information gives AT^ = (0.01/0.7)Amrf and AF^ = 
(0.25/20)Ams, so in both cases AF ~ O(10~^)Am. Equation ^ then implies, in 
general, that |Fi2| = 0(10-2)171^121, leaving 

Am = 2|Mi2| and 

AF = 2Re{Mur*u)/\Mi2\. (23) 
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Hence, the final formulae are as follows: 

Arrid = ^mBjl,BB,ilBmlf2{ml/M'^)\V:M\ 

Ams = ^mBjlBB,VBmU2{m^t/M^)\V:M' (24) 

and 

Am^ _ mBjlBB, 
^rud ^bJI^Bb^ 

Assuming unitarity of the CKM matrix, \Vtb\ ~ 1- Since the top quark mass 
has now been measured, the only unknowns in these formulae are the magnitudes 
of the CKM elements Vtd and Vts , the decay constants /s (which in principle can 
be measured) and the bag factors, Bb-, which are theoretical constructs. 

Estimating and Bb for the B^ and i?^, then, is essential for extraction of 
the CKM elements. Since these mesons differ only by their light quarks, it is 
reasonable to expect that model dependences in predictions of, for instance, Jb^ 
and /^^ are significantly reduced in the ratio fBa/fBs- These issues are addressed 
the next section. 



Vtd 



(25) 



2.3 The B Decay Constant and Bag Parameter 

The B decay constant /b is necessary to relate the B oscillation frequency to the 
CKM matrix element Vtd, as described in Section ^. Prospects for measuring this 
decay constant directly are rather bleak. In principle, it can be measured just as 
/tt is measured, using leptonic decays. The partial width of the decay B Ivi is 

i2 IT/ |2 



T{B^^l+vi) = ^' /ImBmf 1-^ , (26) 

ovr \ mi 



where I = e, fi, r. The decays into e and ^ are helicity suppressed, leaving the 
experimentally challenging B — > tVt- as the dominant leptonic decay. The branch- 
ing ratio into this mode, however, is expected to be very small. Using the range 
for Vud of 0.002-0.005 @ and Jb^ = 200 MeV (see below) results in a branching 
ratio between 1.8 x 10~^ and 3 x 10~^. 

ALEPH 0, DELPHI pl and L3 |l|] have placed 90% CL upper limits on 
this branching ratio, of 1.8 x 10~^, 1.1 x 10^"^ and 5.7 x 10~^ respectively. Since 
the LEP data is played out, no further improvement is likely. CLEO searches 
for a S — > TUr signal by reconstructing the mass of the B meson, using the beam 
constraint to account for the missing neutrinos. The background to this analysis 
is very high - the signal is —9 it 36 events. The CLEO group reports a 90% CL 
limit of 2.2 x lO'^ based on a sample of 2.2 x 10^ B decays. 

The most promising method to extract /b from experimental data had seemed 
to be via the B^ — B^ lifetime difference. It was recently estimated [^] that the 
lifetime ratio is 

with an extra 0(15%) error added to account for terms of order mf that are not 
included in the calculations. Thus, a measurement of the lifetime ratio to 1% 
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corresponds to a measurement of /b to ~ 20 MeV, with a 30 MeV theoretical 
error. The current value of this ratio is 1.04 it 0.04 M. CDF has measured it as 



1.06±0.07±0.01 in exclusive B J/^K decays g§ and as 1.11±0.06±0.03 in 
semileptonic decays |23]. The J/^ sample in CDF's upcoming Run II is expected 
to be more than 20 times larger than the current one, with a similar increase in 
the inclusive sample. Along with modest improvements in the systematic errors, 
this leads to an expected precision on this ratio of « 1%. The BaBar group 
expects to measure this ratio to close to 1% as well |24]. Recently, however, 
our ability to predict lifetimes and lifetime ratios has been called into question 
because none of the models account for the short lifetime. This failure has led 
some authors [25| to suggest that some of the approximations used in deriving 
Equation ^ are not valid - specifically there are large nonperturbative effects 
that should not be ignored. Until the situation is clarified, extracting physical 
quantities from heavy quark lifetimes is suspect. 

This leaves us with two ways to obtain Jb'- lattice QCD and QCD sum rules. 

In the past year or two, lattice calculations of Jb have matured considerably. 
Our understanding of lattice discretation effects has grown, and the advent of 
nonrelativistic actions have provided a means to greatly reduce the uncertainty 
introduced by the relatively large lattice spacing. 

One of the main problems with calculating Jb on the lattice is that the h 
quark is very heavy, and the wavelength of the b quark is small compared to the 
lattice spacings currently computationally feasible. Typical lattice spacings, a, 
available today are 0(0.1 — 0.2) fm, so one has mqa ~ 25a > 1 (where a is in 
fm). Thus, in extrapolating to the continuum limit, mqa = 0, one has a long 
way to extrapolate; further, one must contend with lattice-size artifacts which 
are 0{mQa), ©[(mga)^] and so on, and each may be sizeable. 

The original Wilson lattice action was developed for light quarks, for which 
ma <^ 1 holds. This discretization procedure introduces lattice spacing artifacts 
of 0{a), 0(0^) and so on, which makes the extrapolation to the continuum, 
a = 0, a major source of systematic uncertainty. The Sheikholeslami-Wohlert 
(SW or clover) action [^] was introduced as an improvement on the Wilson 
action for light quarks; it removes the leading 0{a) lattice effects, resulting in 
a much-improved behaviour. Lattice calculations continue to use both of these 
actions for the light quarks in the B and D meson systems, along with a variety 
of approaches in handling the heavy quark. 

In placing a heavy-light system on the lattice, different techniques were devel- 
oped for handling the heavy quark. Taking the straight Wilson action is prob- 
lematic because mqa > 1. One common method of dealing with the b quark, 
known as the static approximation, is to set its mass to be infinite. In the first 
six or so years of lattice calculations this method was used extensively, and it 
produced results for Jb which varied over an extremely wide range. By 1994, 
these variations were understood to be due to large contamination from excited 
B* states in the lattice version of the B meson. Once this was corrected for, 
the results from simulations which used static b quarks stabilized, with values 
comparable to those from other methods. 

In the meantime, propagating b quarks were introduced by performing the 
lattice calculations with finite quark masses, but closer to nic, where mqa < 1. 
Typically, several calculations with different masses were performed, and then 
used to extrapolate up to mq = m^ using a linear extrapolation in 1/mq. Often 
the static limit was included to help anchor this process. In these methods. 
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systematic errors are introduced both by finite lattice spacing artifacts and the 
uncertainties in the mass extrapolation. 

More recently the heavy 6-quark problem has been addressed by using an ef- 
fective nonrelativistic action for the b quark. In the rest frame of a i? meson, the 
light valence quark has a momentum ©(Aqcd)- By momentum conservation, the 
heavy b quark must also have momentum of this order. If the b quark is nonrel- 
ativistic, p = tuqv = Aqcd, and so the velocity v AqcD/^iQ ~ 0.05, justifying 
the nonrelativistic approximation. Thus, one can effectively expand quantities in 
powers of l/mq and expect reasonable behaviour. Further, in a nonrelativistic 
effective theory, the rest mass of the b quark does not appear, so we sidestep the 
problem of having a large mqa. Since the first use of this approach in 1994 [^], 
which included the ©(l/mg) operators, various groups have continued to add 
more and more of the 0(1/771^) terms. 

In such nonrelativistic QCD (NRQCD) calculations, l/mq effects and lattice 
spacing a effects end up intertwined, so one of the challenges of such work is 
to effectively evaluate systematics due to the extrapolations. The uncertainties 
are much reduced if the next-order l/rriQ terms are included. In early 1998 the 
MILC collaboration presented results with a comprehensive treatment of both 
these effects, with improvements through ©[(Aqcd/^iq)^] and 0{a) |4^ . 

A third technique, developed by the Fermilab group (KLM) (5^], allows the use 
of either Wilson or clover actions, with suitable normalisation and nonrelativistic 
interpretation, for the quarks at any mqa. Simulations based on this heavy-quark 
action are run directly at m = m?,, so no mass extrapolation is required. 

Until recently, all calculations of Jb have been performed within the quenched 
approximation, i.e. no virtual quark loops were simulated. Estimating the size of 
the shift in fs due to this omission has been one of the great outstanding issues 
for the lattice community. There are now three groups reporting partially un- 
quenched simulations, with two flavors of light fermions: the MILC collaboration 
pO| , |55[ |, Collins et al [51|, and CP-PACS |52, |^]. The results are summarized 
in Table ^. All three groups found that unquenching the lattice raised all the B 
and D decay constants by 15-20%, although the results are still preliminary. 



Table 3: Results for the decay constant Jb^ from several lattice simulations with 



two dynamical sea quarks. 


Reference 


fB, (MeV) 


MILC 99 [50| 


194 ± 22^^^ 


CP-PACS 99 (5 


\ |5| 210i?6 ± 20 


Collins 99 |l| 


186 ± 25^^° 


Average 


200 ± 30 



Figure ^ shows lattice results for fs^ over the 10-year history of the field. 
The results are grouped such that similar methods of dealing with the b quark 
(static, mass extrapolation, KLM, NRQCD) are shown together for easier com- 
parison. The points are grouped into three rough time periods: pre-1994 years, 
the early years; 1994-1997, during which NRQCD techniques were developed to 
include higher-order terms; and post-1997, when fully 0(1/771,^), 0{a) improved 
simulations were implemented. 
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Figure 4: Summary of lattice calculations of the decay constant /s^. The re- 
sults are grouped by the treatment of the heavy quark: down-pointing triangles, 
mass extrapolation from m < mc, up-pointing triangles, static calculations; open 
crosses, mass extraploation from m > ra^ to rrih] circles, NRQCD; squares, the 
KLM action at rrif,. The rightmost three results, marked with a star, are those 
with two dynamical fermions. 



Table ^ shows averages of fs^ , foa > and fs^ I fsd > using data from the collabo- 
rations shown. Results with full estimations of systematics were included in the 
fit. The average result, /b^ = 166 ± 8 MeV, does not include a systematic due 
to quenching. This is much lower than the average obtained from the A'^^ = 2 
calculations, which is 200 it 30 MeV (allowing a large systematic as these are still 
preliminary) . 

Unfortunately, checking the lattice and sum rule results against experiment is 
difficult, because experimental data are in short supply as noted earlier. However, 
there is one decay constant prediction which can be correlated with data, namely 
/Ds- The branching ratio of Ds — > tu^- is fairly large, 7 ± 4% Q, owing to the 
relative lack of hadronic decay modes in contrast to both B and D decays, which 
have ample phase space available for a multitude of hadronic modes. There are 
currently eight measurements ||5^ of the branching ratio Bv{Ds — > tu-t), which, 
along with Vcs, allow us to extract Jd^ as per Equation |2^. Figure ^ summarizes 
these measurements. 
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Table 4: Lattice results from simulations selected as described in the text. A 
common systematic error due to quenching is not included. 
Reference /b, (MeV) /p (MeV) /bJ/b, 

MILC 98 i6i 157lfl 2101^;^ 1.11 ±0.05 



JLQCD 98 gH 173 ± 13 224 ± 19 1.18±0.06|3| 
APE 98 ill 1791^^ 231 ± 13 1.14 ± 0.03 



GLOK 98 p|] 147 ± 19 NA 1.20 ± 0.04 

FNAL 98 Jp] 164 ± 16 213 ± 18 1.13 ± 0.05 

APE 97 180 ± 32 237 ± 16 1.14 ± 0.08 

Average 166 ± 8 226 ± 8 1.15 ± 0.02 ± 0.06^ 

^ Following a recent review [pi] a common systematic of 0.06 is assigned. 



The ARGUS result is model-dependent and omitted from the average. Results 
from WA75 and E653 have been adjusted to take into account updated Dg and 
branching ratios. Note that the errors have actually increased as the Ds — > (pir 
branching ratio measurement was downgraded by the PDG from having an error 
of 0.4% in 1994 to 0.9% in 1996. The measurements were combined, taking into 
account correlations in the branching ratios used to derive ■ The result is 

f^P = 273 ±28 MeV. (28) 

We obtain a lattice average of fn^ = 226 ± 8 MeV, not including quenching 
errors, from the same group of collaborations used to derive Jb^, as shown in 
Table ^. There is reasonable agreement between the measurement and the lattice 
result, though the agreement would be better using the unquenched lattices. 
More precise experimental data and a better understanding of the unquenched 
lattices would be of great help. 

Since is likely to be the only experimental input for many years, an in- 
teresting exersize is to bootstrap fs^ and fs^ from it. The idea is that many of 
the largest systematic errors in the lattice calculation of the heavy-light decay 
constants either cancel, or are largely ameliorated, when one takes a ratio. For 
example, unquenching the lattice is expected to have a similar effect on Jb^ and 
) as should lattice size corrections. Several groups gs], ^ are now calculating 
the ratio fB^/fOs directly from the same simulation. Multiplying the average 
value of fBa/ fos = 0.78]t4 by the experimental value for results in: 



^lat 

Je 

fh 



fB, = ^/S:P = 213±22(expt.)l}?(theor.) (29) 



The other technique used to calculate Jb is based on QCD sum rules. (For an 
excellent summary of this method, see Reference [^^.) The sum rules are derived 
from the two-point correlation function: 

n(g2) = i J d^xe''^'' < 0|rg(x)i756(x), 6(0)^75^(0)10 > . 

Inserting a complete set of states with S-meson quantum numbers in the time- 
ordered product gives 

2^ _ < 0\qi^5b\B >< B\bi-f5q\0 > ^ < Olqi'y^blh >< hlbi-y^qlO > 

11(9 j — 2 9 / J 2 9 ■ 

— ^ ml — q-^ 
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ARGUS 92 
WA75 93 K 
BES 95 
E653 96 
L3 97 
DELPHI 97 
CLEO 98 
ALEPH 98 
Average 



-4^ 



267 + 28 MeV(Model) 
209 ± 41 ± 48 MeV 
430±150±40MeV 
200 ± 36 ± 32 MeV 
309 ± 58 ± 47 MeV 
330 ± 95 MeV 
280 ± 34 ± 34 MeV 
284 ±18 ±60 MeV 
273 ± 28 MeV 



Gim98 
BBS98 1^ 
JLQCD98 |l « l| 
JLQCD95 I 
UKQCD96 
BS96 
ELC92 
BDHS88 
Average 
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fosC^eV) 



1 



1.5 



Figure 5: Left: Average of experimental measurements. The model- 
dependent ARGUS result is not included in the average. Right: Average of 
lattice calculations of B^^'^. These have been adjusted so that the NLO correc- 
tions are handled in a uniform manner, and include a common systematic of 0.07 
due to uncertainties in the renormalization constant. 



Using nib < 0\qij5b\B >= m\fB, we see that the first term contains 

The two-point correlation function can also be expanded using an operator 
product expansion (OPE): 

Ii{q^) = ^C,(Q2,/x)Orf(/x), (30) 
d 

where are the (calculable) Wilson coefficients of the operators of various di- 
mension d which appear in the expansion. 

Equating these expressions for n(q^), along with experimental input on the 
expectation values of the operators Od (e.g. from bottomonium and charmonium 
masses), allows the extraction of fs- 

The OPE expansion of n(g^) is rather sensitive to the value of the 6-quark 
mass chosen, which ends up being the dominant systematic effect. For example, 



taking m;, = 4.7 ± 0.1 GeV results in /s = 180 ± 30 MeV |7|. 

Given that unquenching the lattice seems to shift /^^ by 15-20% and that 
several groups confirm this result, we choose to use the preliminary results with 
Nf = 2 rather than perpetuate what seems to be a low /^^ value. The quenched 
results on Jb^ measured in a multitude of ways seem to have stabilized, which 
suggests that the unquenched simulations will soon have errors under much better 
control. 

Combining the three values: 

fsa = 200 ± 30, (Lattice, Nf = 2) 
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fs^ = 180 ± 30, (QCD sum rules) 

fB, = 213±22(expt.)t}?(theor.) {fil/fDlxfoD 
gives our estimate: /b^ = 200 ± 16 MeV. 

2.3.1 THE BAG PARAMETER 

The second quantity standing between measurements of Am^ and the clean ex- 
traction of Vtd is the so-called bag parameter Bb- Bb represents a "fudge-factor" 
introduced because we collapsed the sum over all intermediate states in Equa- 



tion 20 to just the vacuum intermediate state. If the vacuum is the only state 
which ends up contributing, it saturates the sum, and Bb = 1- 

As such, Bb IS a, completely theoretical construct and cannot be measured ex- 
perimentally (only the product /b^b has physical meaning). Bb is also amenable 
to study on the lattice, and many estimates of its value have been made by the 
same groups which calculate /b- 

To calculate Bb, the ratio 

<B\dLL\B> <B\dLL\B> 



< B\6ll\B >vac 8/3/1^1 

is evaluated using standard lattice techniques (where Vac refers to the vacuum 
saturation approximation). As with Jb, various groups have used static b quarks, 
propagating b quarks with mass extrapolation, and NRQCD actions. Bb depends 
on the choice of scale fi, and groups quote their results at scales ranging from 2 
GeV to 5 GeV. The lattice and continuum operators must then be matched at 
this scale. The physical quantity of interest ~ the NLO, renormalization-group- 
invariant quantity, is 

52LO(/.) = a.(^)-2/^«(l + ^J)S^(M), (31) 
where 

J ^ Toft. _ J7]_ = 4 

^ 2/3o' '0 
a 11 2nf „ 4nf 

Po = 11 "V^' 71 = -^+9 

= 51-ifi 

Throughout the literature, a number of different choices are present, e.g. nj = 

0,4,5, Ag^o = 200,239 MeV, A§J,j3 = 183,200 MeV, mf, = 4.33,5.0 GeV, as 
well as scaling with J = 0. 

Following Reference |59| , we have tried to rescale the results given in a consis- 



tent manner; namely, n/ = 4, Aq^j-, = 200 MeV, and = 5 GeV are used to scale 
CKs up to the matching point at 5 GeV, with Og given by the NLO expression |^] . 
We match between nj = 4 and nj = 5 at the b mass of 5 GeV, and set Uf = 5 
in the exponent of q^, —2/Pq. This results in a scaling of about 1.6 between the 
raw BB{mi)) and B^^'^{m},). 

The results on Bb^ from calculations using static b quarks were omitted from 
recent averages due to an extremely large variation in the results when the NLO 
corrections were included in three supposedly equivalent ways pi, lei. Recently, 



several clarifications to this process were made by Gimenez & Reyes |5^. First, 
the calculation of the lattice operator Oll = ^7^(1 — 75)'7&7^(1 — 75)9 for SW 
(clover) actions used previously was incorrect, except for those groups using 
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Wilson actions for the light quarks. Second, the importance of using tadpole- 
improved light quarks was shown. Third, some inconsistencies in the techniques 
used to match operators to the continuum were pointed out. They have rean- 
alyzed their data along with data from the UKQCD group |65|, which also 
used SW light quarks. The results are now much more stable under variations 
in the matching scheme, as well as being in good agreement with the data from 
simulations that used propagating heavy quarks |^, 64, |66[ . 

Figure ^ shows the average for B^^^ from the analyses in References [28, 33, 



62 



65, 



which is 



1.37 lb 0.08. A common systematic error 



of 0.07, due to uncertainties in the renormalization, has been included. 

The ratio Bb^/Bb^ is expected to be very close to one, and has been estimated 
as BbJBb, = 1.01 ± 0.01 H" 



3 ANATOMY OF A TIME-DEPENDENT OSCILLATION MEA- 
SUREMENT 

We now turn to the experimental side of explicitly measuring the time depen- 
dence of the BB flavor oscillations. This section summarizes the basic analysis 
elements and techniques that have been used to produce the results presented in 
the following two sections. 

The basic steps of any time-dependent mixing measurements are evidently as 
follows: 

1. Find the interaction point at which the B mesons were produced (primary 
vertex) and decay point of one of the mesons (secondary vertex), and deter- 
mine the decay distance L. 

2. Measure the (presumed) B meson's momentum to turn the decay distance 
L into the proper time lived, t. 

3. Tag the flavor of the (presumed) B at its production point, and determine 
the probability the tag was correct, P^{t). 

4. Tag the flavor of the B at its decay point, and determine the probability 
the tag was correct, PP{t). 

The advent of high-precision silicon tracking systems allows the measurement 
of the decay distances of B mesons, which are typically of order "f(3cTB ~ 7/5 x 
470/im. With boost factors ranging from /37 ~ 3 at CDF to /37 ~ 7 at LEP and 
SLD, decay lengths in the millimeter range are typical, which are easily discerned 
by modern detectors. For detectors with three-dimensional precision tracking (i.e. 
double-sided silicon or pixel detectors) t = L/^j3c = LMb/pb- Experiments with 
single-sided silicon strip detectors measure the flight distance projected onto the 
plane perpendicular to the beamline (xy, or r(j) plane), L^y, and convert to L 
using the estimated polar angle of the meson, t = Lxy/i'j (3c sin 6) = L^yMB/pt- 

The uncertainty on the decay time of the meson in its rest frame is 

^(t) = + (32) 

The resolution on L (or L^y) is typically C'(300//m) at LEP and ©(lOO/um) 
at CDF, so when divided by 7 does not contribute significantly to the error 
on ct compared to the momentum-resolution term 17(7/3), which grows linearly 
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Figure 6: Schematic of a typical B event. 



with the decay length L. This is a critical point for measuring B^ oscillations, 
which have a very rapid oscillation rate. Once the smearing induced by the 
momentum resolution reaches half the oscillation period, any hint of an oscillation 
signal is washed out. With the momentum resolution achievable from partially 
reconstructed decays, only the first few oscillations of the B^ will be visible, 
making frequency measurement very difficult. Fully reconstructed B^ decays, of 
course, have optimal momentum resolution, but current experiments have very 
small data samples. 

The variety of mixing measurements comes mostly from the wide range of 
techniques used to identify the 6-quark flavor at its production and decay point. 
These techniques are breifly described below. 

By far the dominant production mode of 6-type quarks is bb pair production, 
from e^e~ annihilation at LEP and SLD, and from quark- antiquark annihilation 
and gluon-gluon fusion at CDF. At collider experiments they are generally pro- 
duced roughly back-to-back, so the event can be divided into two hemispheres 
(defined by plane perpendicular to a vector a passing through the event origin) 
that each contain one of the b quarks. The analyses are insensitive to slight 
changes in the direction of d, and several choices for a are used, including: the 
thrust or sphericity axis of the event; the direction between the primary and 
secondary vertices; or the direction of a jet associated with the secondary ver- 
tex, defined by an appropriate clustering algorithm. The hemisphere with the 
secondary vertex is called the vertex side, and the other is called the opposite 
side. (This distinction becomes arbitrary if vertices are reconstructed in both the 
b and b hemispheres.) 
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Figure ^ is a sketch of an event, showing the b and 6 quarks originating from 
the primary vertex, a secondary vertex found on the b (vertex) side, and some 
possible flavor tags on the opposite side (see below). The hemispheres would be 
formed by a plane passing through the primary vertex roughly perpendicular to 
the b-b axis. 

3.1 Decay-Point Flavor Tagging 

The 6-quark flavor at decay is typically identified by one of two methods, either 
by the charge of the lepton(e or ^u) from a semileptonic decay of the meson, or 
by the type of D meson present in a partial reconstruction of the decay. 



Table 5: Sources of leptons from b quark decay (correct flavor tags of the b quark 



are denoted r=right, 


incorrect tags are denoted w=wrong). 


Decay chain 


Lepton Source 


Correct Tag? (r/w) 


b^W- -^l-, 


Direct 


r 


b^c^W+ ^1+ 


Cascade 


w 


b^W- ^c^l~ 


Right-sign cascade 


r 


b^W- ^T- ^l- 


Right-sign cascade 


r 


b J/i)X /± 


Both-sign cascade 


r/w 



A lepton from a semileptonic b decay can have the sources given in Table |. The 
main contributions come from the first three components. The relative rates of 
leptons from direct and wrong-sign sequential decays are approximately equal |p 
The right-sign sequentials 6 — > c — > /~ contribute approximately the same amount 
as the wrong-sign Q, so the lepton tags the 6-quark correctly approximately 2/3 
of the time. However, with very simple kinematic cuts, the contribution from the 
sequential decays can be greatly reduced. In particular, the momentum spectrum 
of the sequential leptons is softer than that of the direct leptons, as they arise 
further down the decay chain. 

The second common method of tagging the decay-point b charge is by either 
partially or fully reconstructing the charmed meson from the B decay. In par- 
ticular, the decays W D*+X, {D*+ D^-n+) and W D+X allow the 
separation of b from b quarks, (vr,,, denotes the soft pion from the D* decay.) The 
sign of the 6-quark charge is the same as the sign of the kaon charge in the D 
meson decay, and it is opposite that of the tt^: from the D*~^ decay. 

Figure ^ shows the decay modes of the B'^ and B~ assuming final-state in- 
teraction effects are small (i.e. factorization approximately holds). Soft gluons 
running around the diagrams are not shown. In the semileptonic modes (Fig- 
ure @a,&), note that the B~ produces a D'^ , D*^ , D**^ and that the D*'^ decays 
only to D^, not D'^. Thus the only contamination of B~ in the D*~^ sample 
is through decays involving D** . In contrast, the B^ decay can produce a 
directly, or through a D*~^ with a missing vr*^ and a D*~^ detected in the D^nf 
configuration. The decays involving D** will end up scattered around the various 
modes, as in the B~ case. 

Hadronic decays show the same basic behavior, with the addition of the inter- 
nal, color-suppressed diagrams (Figure 0e,/). Note that the type of charm decay 




Figure 7: Diagrams for semileptonic and hadronic decay of B'^ and B mesons 
(soft gluons are suppressed). 



product for these diagrams is the same for the B~ decays. The B^, however, 
would tend to produce D^,D*^, identically to the B~ . This breaks the anal- 
ogy between the semileptonic and hadronic modes, resulting in a slightly different 
fraction of D*~^ coming from in the two cases, but the effect is not large. 

Diagrams like those in Figure and can also have the external W decaying 
into a cs pair, producing a charm of the same charge as the b quark. Fortunately, 
the branching ratio of B ^ L>(*)I)7^*^ is only 0(5%). 

Evidently, requiring a reconstructed D*~^ (with or without an associated lep- 
ton) in the final state not only flavor-tags the b quark at the decay point but 
also greatly enriches the sample's B^ fraction. Although no measurement of the 
fraction of D*~^ that comes from B^ (as opposed to B~) is available for hadronic 
decays, it has been measured to be, e.g., 84 it 9% by OPAL in the semileptonic 
mode. 

The advantages to the D*~^ tag are its excellent purity and its ability to produce 
event samples with very high B^ content. Its obvious drawback is its efficiency, 
which is rather low owing to the exclusive final state. 
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A third method of tagging the 6-quark flavor at the decay point is via the 
charge of the kaon(s) in the decay chain. The decay of the b quark proceeds 
6 — > c — > s. Hence a K~ tags a b quark and a tags a b. Using inclusive kaons 
will greatly increase the efficiency of this tagging method over the D* method, 
but of course there are other sources of kaons besides the Z?-meson decay, which 
will dilute the tag's effectiveness. In addition, the dramatic increase in the 
content of the sample is lost. 



Finally, a method based on a charge dipole is used by SLD (see Section 4.6). 



3.2 Production-Point Flavor Tagging 

The flavor tagging of the b quark is more problematic at its production point 
than at its decay point. Broadly speaking, there are two strategies: same-side 
techniques, which use information from the vertex hemisphere to determine the 
initial flavor of its b quark, and opposite-side techniques, which measure the 
charge of the b quark in the opposite hemisphere (at either its decay or production 
point) and infer the vertex-side b charge from that. 

Opposite-side tags come in two types: lepton tags and jet charge. Lepton 
tagging is essentially the same as described above for decay-point flavor tagging, 
although the selection criteria on the lepton may differ for the two sides. This 
type of tagging has high purity but suffers from a low efficiency because the 
semileptonic branching fraction of B mesons is only 20%; in addition, detection 
and selection criteria are imposed. 

The jet-charge method relies on a weak correlation between the momentum 
and charge of tracks in a bb event. For example, at SLD or LEP, b and b quarks 
are generally produced back-to-back each with momentum of 45 GeV. Thus the 
b quark hemisphere starts off with a net —1/3 charge. As the fragmentation 
proceeds, the b quark gradually loses its momentum as it radiates off gluons; 
however, the —1/3 charge continues to be associated with a higher-momentum 
particle, whereas the charge required to keep the hemisphere charge integral 
is in lower-momentum particles. Experimentally, a momentum-weighted charge 
distribution is used to form a variable to discriminate between b and b quarks. 



Section 4.3 provides details on the algorithms used to measure the jet charge. 

Two types of same-side tagging (SST) are commonly used. One method relies 
on an even more specific charge correlation than the jet charge [^]. Figure ^ 
shows a schematic of a 6-quark fragmentation chain. If the first quark-antiquark 
pair radiated is a uu, the b becomes a B^ , and the first charged pion in the 
fragmentation chain is a vr". If the first qq pair is a dd, the meson is a B^ and 
the first charged pion is a vr"^. Hence for a B^ / B^ sample, a correlated pion vr^ 
tags a B^ and 7r~ tags a B'^. The correspondence is opposite for charged B^. 
This tagging technique, then, requires a rather clean sample of B^/B^ events, 
which rules out applying it to purely leptonically tagged events. The difficultly 
lies in selecting the correct pion from all the fragmentation tracks. In the case of 
a fragmentation into a B^ meson, the first fragmentation track would be a K~ . 

The second same-side method is unique to SLD, and relies on the large lon- 
gitudinal polarization that the SLC was able to achieve. The differential cross 
section for producing a b (as opposed to b) quark in a polarized e'^e~ interaction 
is 

abicOsO) = {I- AePe) {l+COs'^e) + 2 Ab{Ae-Pe) cos 0, 

Gi COS (7 
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Figure 8: Schematic of and B~^ fragmentation chains 



where 9 is the angle between the b quark and the incoming electron direction. Pg 
is the longitudinal polarization of the electron beam (Pg > for positive helicity, 
right-handed polarization), and Af = 2vfaf/{v'j +a^) where af,Vf are the axial 
and vector coupling of the fermion / to the Z^. Thus the probability a ^'-type 
quark produced at angle 9 is a b, rather than b is simply 



cjfe(cos 9) 



Ab{Ae- Pe) cos ( 



(Tb(cOS 9) + o"^(cos ( 



1 

2 (1 -^ePe)(l +C0S2( 



= Ppol(cOS( 



(Some papers have the sign of Pg swapped and/or measure 9 from the positron 
direction.) The polarization of the electron beam is randomly chosen between 
the positive and negative helicity states on a pulse-by-pulse basis, and the signed 
polarization is recorded in the data stream. The standard model values for the 
asymmetries are A^, = 0.155 and Af, = 0.94 Q. Note that the discriminat- 
ing power that survives at Pg = because of the residual forward-backward b 
asymmetry has a maximum amplitude of AgAf, ~ 0.1, whereas for large Pg this 
amplitude is approximately PeAi, ~ Pg- Hence SLD can tag the flavor of the b 
quark at production simply by measuring the direction of the P-meson system. 
The probability that the sign assigned is correct, for a given angle, is simply the 
-fpoi(cos6') above. It is the large (63 and 77%) polarizations achieved at the SLC 
that make this flavor tag effective, with a mean right-tag probability < Ppoi > of 
62% and 76% respectively. Figure 10 shows a typical distribution of cos 9, along 
with the underlying predictions for b and b quarks given by the Monte Carlo 
simulations. 

Finally, many experiments combine multiple tags, both same-side and opposite- 
side, to maximize the probability of tagging correctly. For example, the SLD 
group uses the event-by-event Ppoi in all their analyses, along with a second flavor- 
tagging technique (e.g. jet charge). The two tags, Pi and Ppoi, are combined into 
the final right-tag probability 



P 



PP 



comb 



l^pol 



^'li^pol + (l-^'l)(l-i^i 



poU 



(33) 
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3.3 Fitting Procedure 

The signature for S-meson mixing is that the sign of the b quark inferred at the 
decay point, D, is opposite to that inferred at production, P. In the fohowing, the 
notation DP+ refers to events in which the experimental flavor-tag signs for the 
vertex-side b quark at production and decay are the same, (i.e. the event fooks 
unmixed) and DP— is used when they are opposite (i.e. the event looks mixed). 
For example, in an analysis that uses leptons to tag both the production and decay 
charges, DP+ events are those with opposite-sign leptons. The event sample is 
then split into two components, with differing time dependences Nop+it) and 

NDP-it)- 

The basic variable to be measured, then, is the time-dependence of the decay 
and production flavor tags' charge correlations, qdt) = sign{DP)(t). The mean 
value of this variable, binned in time t, is the standard measured asymmetry: 

Let us first consider this time-dependent charge correlation on a sample of pure 
/ B^. The numbers of events in the DP+ and DP— samples at a given time 

are 

NDP+{t) = N(Pdt)[PfP,'' + {l-Pf){l-Pf)] 



+ Pm{t)[Pf{l - P?) + (1 - P^)P?]) 

NnP-{t) = N^PUmPfPr^' + il-PfX^-Pf)] 
+ P„(t)[P/'(l-P,^) + (l-P/')P,^ 



where is the total number of events in the sample, Pu and Pm are the unmixed 
and mixed probabilities from Equation 14, and P^{PP) is the probability that the 



production (decay) flavor tag got the sign of the b quark correct. The asymmetry 
is 

Am{t) = (2P/'-l)(2n^- 1)^11^^ =PpPz.cosAmt. 

That is, the measured asymmetry Am is reduced from the true asymmetry A{t) = 
cos Ami by the production and decay tag dilutions, Vp = 2P^ — 1 and Vd = 
2P^ - 1. 

Now suppose we have a variable q, calculated event-by-event, on which the 
production tag probability depends. If we break up the sample into bins of a, 
with N{a) events in a given bin, then (setting P^ = 1 to reduce algebraic clutter) 
we have 

NDP+{t,a)-NDP-{t,a) = iV(a)(2P/'(a) - 1) [P„(t) - P^(t)] 



The error on this difference is \jN{a), independent of the dilution P(a), while 
the value of the difference scales with T){a). 

If we perform the time-dependent fit separatly for each a bin, the relative 
errors on the fits will clearly scale as a{a) ~ l/[(2P,.(a) — N{a)\. These (sta- 
tistically) independent measurements are optimally combined using the weighted 
mean. This is equivalent to weighting each event in the sample by its dilution, 
2P^{a) — 1, and fitting once on this weighted sample. A similar result holds for 
the case with both P^ and P^ depending on a set of event-by-event variables a. 
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Thus in a time-binned fit to the asymmetry, the optimal procedure is to weight 
each event's contribution by the estimated dilution for that event. Most of the 
analyses described here use an unbinned likelihood fit, applied to the N£)p^{t) 
and A'^DP-(i) samples simultaneously. In this case, employing the event-by-event 
right-tag probability (a) is optimal. 

The likelihood functions for the DP+ and DP— subsamples, or the binned 
asymmetry, are constructed based on predicted proper time distributions of the 
data sample's components. These typically include unmixed decays, mixed 
B^ decays, i?^, h baryons, unmixed B^, mixed i?^, cascade decays from all the 
above h sources, charm background, and fake or combinatorial background. 

The measured decay time in an event differs from the true decay time because 
of imperfect reconstruction of the decay length and of the B momentum, as shown 



in Equation 32. The probability of reconstructing a decay time t given that the 
true decay time was t' is described by the resolution function R{t,t'). This is 
convoluted with the predicted true distribution, e.g., Pu{t) or an exponential, 
to get the expected measured distribution. The vertexing algorithm might also 
introduce a f-dependent efficiency (actually L-dependent, which roughly trans- 
lates into i-dependent), which is sometimes accounted for implicitly in R, or by 
multiplying the above convolution with an efficiency curve. 

R is typically parameterized as a sum of two or three gaussians, with sepa- 
rate parameterizations over several bins of decay time to allow for the increased 
smearing at longer times. If the vertex efficiency is included, the gaussians can 
differ above and below the mean to fit, e.g., a decrease in the probability to 
reconstruct a vertex near the primary vertex. 

Alternatively, some analyses fit directly to the predicted distribution as a func- 
tion of decay length, in which case the resolution function is parameterized as a 
function of decay length rather than time. 

4 fiO MIXING 

The time-dependence of B^ — B^ oscillations was first observed by ALEPH [^] . 
Since then, a number of analyses have measured the oscillation frequency, as 
summarized in Table ^. This section presents an overview of these measurements 
grouped by the techniques used to flavor tag the initial and final charges of the 
6-quark. Several preliminary results have been included, which are indicated on 



the summary in Figure 12. 



4-1 Detector Overviews 

What follows is a brief overview of the six collaborations' detectors, focussing on 
the elements which are most relevant to mixing analyses. More detailed descrip- 
tions are available for ALEPH fT^], CDF |7l|l, DELPHI @, L3 |7|, OPAL @] 



and SLD |7|. 

All of the detectors operate at colliders and have similar cylindrical geometries 
with some type of end plug to cover most of the solid angle. Five of the experi- 
ments (ALEPH, DELPHI, L3, OPAL at LEP and SLD at SLC) operate at e+e" 
machines running at the Z-pole, y/s = 91 GeV. CDF operates at the Fermilab 
pp collider, which has an energy = 1.8 TeV. 

Time-dependent oscillation measurements clearly require precision tracking de- 
tectors, and it was only after the inclusion of various silicon strip and pixel devices 
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Table 6: Index of the various flavor-tagging combinations used by each experiment 
to measure mixing. See the references indicated for details on individual 
analyses. 



Vertex tae 



JetQ K/dipole 



Prod. Tae 



JetQ 



/ JetQ 



JetQ SST 



JetQ 



JetQ 



ALEPH 
CDF 
DELPHI 
L3 

OPAL 
SLD 



7t 



781 



[31] 
[32] 



[76] 
[75] 

[31] 

mm 



re 



sg] 



85] 

m 



m 



that each experiment was able to perform these analyses. ALEPH, DELPHI and 
OPAL had silicon detectors in place since 1991. The ALEPH detector measures 
both the r(j) and rz views - OPAL(1993), DELPHI(1994) and L3(1994) soon 
added vertex detectors which could do the same. SLD had a CCD pixel device 
available since 1994. 

Throughout the following, d denotes the impact parameter, or distance of clos- 
est approach, of a track to the beam direction. The momentum p of a track 
projected onto the rcj) plane is denoted pt- 

4.1.1 PRIMARY VERTEX RECONSTRUCTION 

The position of the interaction point of an event (primary vertex) is an important 
part of determining the decay length of the b particle. The dimensions of the 
beam cross section (i.e. in the x — y plane) are small for all the accelerators 
used in these analyses: 150/im x lO/zm at LEP, 25/im x 25/im at Fermilab, and 
1/xm X Ifim at the SLC. Each experiment tracks the mean position of the beam 
centroid, giving an accurate estimate of the interaction point. In addition, some 
experiments fit for a vertex event-by-event, combining this mean position with 
tracks which have a high probability of originating there. This helps reduce the 
relatively large uncertainty on the x position of the primary vertex at the LEP 
experiments. 

4.1.2 PARTICLE IDENTIFICATION 

The analyses presented here rely heavily on detecting and tagging B mesons 
through their semileptonic decays to electrons and muons. These are identified 
using standard detector elements such as tracking chambers, electromagnetic and 
hadronic calorimeters, and muon chambers (details can be found in the refer- 
ences). All experiments remove electron candidates from pair conversion. 

The SLD and LEP experiments have data acquisition systems which can record 
virtually all of the relatively low-rate hadronic Z events. CDF, however, operated 
with an interaction rate of ~ 300 kHz and simply could not record every event; 
only those passing a multistage trigger were selected for analysis. This trigger 
could identify those b events that contained a semileptonic decay by detecting 
the lepton. However, the fake background at low momentum is very large, and a 
fairly high lepton pt threshold of 7.5 GeV was required so that the trigger would 
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not saturate the data acquisition bandwidth. Dedicated dilepton triggers, either 
/i^ or e// were also implemented. For dileptons, the momentum cuts were set 
much lower (2.0 GeV for the muon and 5.0 GeV for the electron), since requiring 
two leptons significantly reduces the fake background. The data sets collected 
with these triggers are referred to as high-pj single-lepton and low-pt dilepton, 
respectively. 

In addition to leptons, fragmentation kaons and pions are used by many ex- 
periments to tag the initial 6-quark charge, and kaons from B decay are used by 
SLD to tag the final charge (see Section |). ALEPH, CDF and OPAL identify 
particles based on ionization (dE/dx) deposits in their tracking chambers, while 
DELPHI and SLD use their ring-imaging Cerenkov detectors. 



4-2 Dilepton Analyses 
4.2.1 EVENT SELECTION 

The dilepton analyses use a lepton tag as the decay-point tag, as described in 



Section |3.1| , and an opposite-side lepton tag as the production-flavor tag. 

Like all opposite-side-type analyses, the dilepton analyses divide candidate 
events into two hemispheres using one of the methods discussed in Section ^. 
Events are selected which contain two leptons {I = e, /x) in opposite hemispheres 
and have a secondary vertex associated with at least one of the leptons. The 
details of the vertexing algorithms are given below. The leptons are required to 
pass the cuts on their total momentum (or momentum transverse to the beamline, 
for CDF) listed in Table 0, which also shows the number of events and number 
of secondary vertices found in each analysis. In addition, most analyses reject 
leptons with a below the threshold given, where pp' is the momentum of 
the lepton transverse to the direction of the associated 6 jet (with the lepton 
removed). OPAL rejects leptons based on the output of a neural net, described 
below. 

CDF has three dilepton analyses - two based on the data recorded with the nfi 
and efi triggers and the third based on events passing the single-lepton trigger 
that have a secondary vertex associated with the lepton. From this sample, events 
with another so-called soft lepton opposite the trigger lepton are accepted (the 
algorithm for selecting the lepton is based on the soft lepton tag, or SLT, used 
in the top quark search). 



4.2.2 SECONDARY VERTEX RECONSTRUCTION 

Because the lepton flavor-tagging method identifies only one possible track from 
the B decay, all experiments use an inclusive vertex-finding technique to locate 
the B decay point. 

ALEPH approximates that all of the 5-decay products except the lepton 
come from the tertiary charm vertex (commonly loosely referred to as the D 
vertex). They form a grid of possible charm vertex positions. If a track's impact 
parameter (excluding the lepton) is within 3a of a candidate vertex, it is assigned 
to the vertex. ALEPH calculates the difference between assigning all tracks 
to the primary vertex and allowing some to come from the charm vertex, and 
the point which maximizes this difference is chosen as the charm vertex. A 
charm pseudotrack is formed, which passes through this vertex, and its direction 
is given by the sum of the momentum of the charged particles in the vertex. 
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Table 7: Summary of selection requirements on the momentum and pi of the 
leptons for the dilepton analyses, along with the number of events and recon- 
structed vertices for each experiment. < > is the mean right-tag probability 
for the lepton opposite the vertex hemisphere. 



Experiment 


Lepton requirements (GeV) 


Events/ Vertices < > 


ALEPH 1^ 


p > 2.0(e), 3.0(/i);pf' > L25 


5957/9710 


0.83 


CDF(/x/i)g 


Pt > 3.0, pf > 1.3 


5968/5968 


0.74 


CDF(e/i) H 


Et > 5.0(6), Pt > 2.5{fi), pf > L25(/i) 


10180/11844 


0.74 


CDF(SLT)#8 


p > 7.5(trig),2.0(tag),pfHweight) 


12700/12700 


0.73 


DELPHI 


p > 3.0, pf^ > 1.2 


4778/4778 


0.89 


L3 g 


p > 3.0{e),4.0{fi),pf > 1.0 


1490/1928 


0.84 


L3(ip) g 


same 


2596/- 


N/A 


OPAL If 


p > 2.0, neural net 


5357/8544 


0.85 



This pseudotrack is then vertexed with the lepton to give the position of the 
secondary, B vertex, which is finally projected onto the 6-jet direction to give 
the decay length L. This technique allows a B decay length to be determined 
for every event and is appropriate for samples with high initial b purity, which 
do not need to rely upon the presence of a well-separated secondary vertex to 
increase the b purity. 

On the other hand, at CDF, the presence of a secondary vertex is a powerful 
means of rejecting large non-b backgrounds. CDF uses two techniques, both of 
which involve the projection of the tracks onto the rcj) plane. The dimuon analysis 
[fr^ j exploits a correlation between the impact parameter, d, of a track, and its 
azimuthal (i;^) angle. Tracks coming from a secondary vertex form a line in the 
d<j) plane, whereas tracks from the primary vertex have both small d and no dcf) 
correlation. A cluster of tracks from a secondary vertex is formed if at least three 
tracks with significant separation from the primary vertex [d/ad > 2) forming 
a line in the dip plane can be found. All tracks in this cluster except the muon 
are vertexed to form a presumed charm vertex, and a pseudotrack is formed, as 
above, from the summed momentum of these tracks. This is intersected with 
the muon track to give the h decay vertex position, and hence the flight distance 
projected onto the r4> plane, Lxy. 

The second method is used for the high-pt-SLT ||7^ and e^u analyses First, 
candidate tracks from the 6 jet are selected if they have significant impact pa- 
rameters to the primary vertex (this is not required for the lepton track which is 
presumed to come from the b decay). Combinations of at least three tracks con- 
sistent with coming from the same point are candidate vertices. If no such vertex 
is found, stricter impact parameter and pt cuts are imposed. All tracks are com- 
bined to a common vertex and tracks which contribute too much to the vertex 
are removed. The process is repeated until a good fit is obtained. This technique 
does not attempt to separate a tertiary charm vertex from the b vertex, but in- 
stead combines all tracks into one. The efficiency of these algorithms necessarily 
falls off steeply as the secondary vertex approaches the primary vertex. 

DELPHI |80| uses yet another method, effective on high b purity samples, 
which attempts to identify tracks from the charm vertex by mass rather than 
position. Excluding the lepton, they cluster particles in the jet associated with 
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the lepton using a jet clustering algorithm optimized for this purpose (all tracks 
are given the pion mass). Within a cluster, particles are ordered by decreasing 
values of their pseudorapidity relative to the cluster direction. Tracks with the 
largest pseudorapidity and p > SOOMeV are kept until the mass of the resulting 
system exceeds 2.2 GeV, and any such system with a large angle (> 500 mrad 
relative to the b jet) is discarded. The remaining tracks in the cluster are vertexed 
(in the r<f> plane), and the charm pseudotrack is formed and intersected with the 
lepton to give the rcj) projection of the b decay point. The decay point with 
the largest significance of separation from the primary vertex is chosen. The 
procedure is then repeated once, using the charm vertex tracks as the cluster 
seed and including neutral particles in the b jet, to obtain the b decay length, 

L3 [^] uses a track-based approach similar to CDF. They search for a secondary 
vertex using the lepton and tracks that are not consistent with coming from the 
primary vertex. An acceptable vertex can be formed in approximately 70% of the 
events, with a loss of efficiency when the secondary vertex is near the primary. 
L3 also has a unique analysis that does not attempt to find a secondary vertex 
at all; rather it looks at the charge correlation of the two leptons as a function of 
the impact parameters to the beam axis (IP) of the leptons. 

OPAL 1 82 1 uses a technique in which tracks, ordered by decreasing significance 
of separation from the primary, are intersected (in the r(j) plane) with the lepton 
track to form secondary vertex seeds. Additional tracks that are more consistent 
with the secondary seed than the primary vertex are added, and a candidate 
secondary vertex is chosen based on its position and the number of associated 
tracks. Quality requirements based on the error on the vertex position, the invari- 
ant mass of the tracks in the vertex, etc., reduce the efficiency to approximately 
70%, though in a decay- length-independent manner. 



4.2.3 B MOMENTUM DETERMINATION 

ALEPH estimates the B meson momentum as follows. First, the (vector) sum 
of the momentum of the lepton and charm vertex tracks is formed. Second, a 
fraction of the neutral energy in the hemisphere, determined from Monte Carlo 
studies to be 0.68, is added to this momentum. Finally, the momentum of the 
neutrino is approximately accounted for by adding the missing momentum in 
the lepton hemisphere, estimated as the difference between the sum of all visible 
energy in the hemisphere and the beam momentum. 

CDF begins with the total pf in the vertex cluster, and scales it to pf/lC, 
where /C is a factor that is determined from Monte Carlo and parameterized as 
a function of the obsevered pf and the mass of the tracks in the cluster. 

DELPHI starts with an estimate of the B momentum given by the total energy- 
momentum reconstructed in its hemisphere minus that of all particles not in- 
cluded as part of the B in the vertex algorithm, which already contains contribu- 
tions from neutral energy. A single neutrino's energy and direction is then solved 
for by assuming it is the only missing particle in the event. The visible energy 
and momentum are rescaled by a factor a = 1.13, which gives the optimal resolu- 
tion. The energy Ei, is added to the B momentum if the direction of the neutrino 
momentum is within 400 mrad of the /-charm direction previously determined. 
Finally this estimate is adjusted to to give the correct mean value, based on the 
Monte Carlo prediction. 
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L3 does not estimate the B momentum event-by-event but rather uses a con- 
stant value of 0.85£'beam- 

OPAL uses a very different technique from the other experiments. OPAL first 
calculates the invariant mass of all the objects in the event that are not in the b 
jet, assuming all charged tracks are pions. This mass is corrected slightly based 
on the inverse of the total reconstructed energy (which should be Mz) to improve 
the resolution. The energy of the b jet, including the neutrino, is determined, 
assuming the undergoes a two-body decay to particles of mass 5.3 GeV. Some 
of this energy is removed as follows. For charged tracks, a weight is calculated 
based on: the probability that the track comes from the secondary, rather than 
the primary, vertex; the track's momentum, and the angle of the track to the 
estimated charm-hadron direction. For neutral clusters, the weight is based on 
the angle between the cluster direction and the charm direction alone. The 
weighted sum of charged and neutral energies in the jet is then subtracted from 
the full b jet energy to give the estimated B meson momentum. This technique 
results in an excellent fractional resolution of approximately 12%, with the core 
being closer to 7%. It was found that the resolution in dilepton events, with two 
missing neutrinos, was not degraded significantly. 



4.2.4 FIT METHOD AND RESULTS 




Proper time (ps) Proper time (ps) 

Figure 9: Fraction of like-sign events vs. proper time. Left: dilepton events from 

ALEPH Right: lepton-jet charge events (selected with IQtagI > 2) from OPAL. 
The curves represent fits to the data. 



The analysis variable for dilepton events is the product of the lepton signs 
Qc{t) = q"^^ ■qj^^{t). Opposite-sign leptons make up the "unmixed" DP+ sample, 
and same-sign leptons the DP— sample. The proper time dependence of these 
samples is predicted using the resolution-smeared, efficiency- weighted distribu- 
tions of the sample components, weighted by their relative abundance, combined 
with the probability of getting the correct charge assignments (as described in 
Section ^^ ). If secondary vertices are found for both leptons, the joint probabil- 
ity of finding decays at their respective decay times is constructed. The CDF //// 
and efi analyses use these predictions in a binned least-squares fit, while the rest 
of the analyses use unbinned likelihood fits. 

The b purity of these samples is fairly high, ranging from 94% for DELPHI, 
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to 98% for ALEPH, with the CDF samples typicahy at 80%. The remaining 
events are either cc, or fakes, cc events represent between 0.5% of the sample 
(for ALEPH, L3 and the low-pi CDF samples) and 8% (for CDF). The fake 
background, including fake leptons from both bb and non-6 sources, constitutes a 
few percent - up to 15% at CDF. Approximately 38% of the leptons in bb events 
are estimated to come from mesons. 

Many of the following analyses use the opposite-side lepton as a production 
flavor tag, so a few words concerning its effectiveness are in order. The lepton 
opposite the vertex side can mistag the b in that hemisphere if the the b hadronizes 
into a B^ or B'^ and mixes, if the lepton comes from a cascade b ^ c ^ I decay, 
or if the lepton is a fake. Denoting the fractions of the various B hadrons by 
their light quarks as fd,fu,fs, and /a^ for all the baryons, the contribution to 
the mistag rate from mixed B^ and B^ mesons is x = fdXd + fsXs (ingnoring for 
the moment those events with vertices on both sides). Since Xd = 0.172 it 0.010 
05 Xs ~ 0-5 (see Section ^), then x = 0.118 ± 0.006 0, with contributions from 
fd ~ 40% and fs ~ 11%. Thus, even though B^ mesons only represent 11% of 
the B mesons produced, they contribute substantially to the overall mistag rate 
due to their fast oscillation rate. 

The other significant source of mistags are the wrong-sign cascade decays. Ta- 
ble |5| lists the sources of leptons which do not come directly from the b quark. 
Among these, the first two make up the bulk of the events. The right-sign prob- 
ability, then, is Pr = l — x~ fb'c^ where /^^ is the fraction of wrong-sign cascade 
decays from the mix of b hadron types in the hemisphere. In the absence of any 
cuts on the lepton, this fraction is comparable to the contribution from direct b 
decays. To reduce this fraction, all experiments try to remove, or assign small 
weights to, events with a large probability of having a lepton from a cascade 
decay. Each experiment uses the momentum of the lepton referenced to the di- 
rection of the jet (with the lepton removed), pl'^^. Leptons from direct b decay 
can be produced with significant mometum perpendicular to the b flight direction 
due to the large 6-quark mass, whereas those from the charm decay typically have 
less transverse momentum. For the same reason, cascade leptons also tend to be 
less well isolated from the core of the jet, and because they come from further 
down the decay chain, their mean momentum is lower. OPAL feeds these three 
variables into a neural net analysis, which is estimated to be 40% more effective 
at rejecting cascade decays than simple cuts on p and These cuts reduce the 
right-sign cascade component, b ^ c ^ I, even more, since these events are of the 
B DD type, and the momentum available for the lepton is even smaller than 
in other cascade decays. These end up accounting for 12 — 25% of the cascade 
decays. 

The SLT analysis of CDF, rather than simply cutting on pp', weights events 
according to the probability that the flavor tag lepton identifies the b charge 
correctly (parameterized as a function of pp^). The other analyses use the mean 
probability for each event. 

The sample composition on the vertex side need not be the same as on the 
opposite side. For example, in the CDF SLT analysis the vertex-side lepton 
has pt > 7.5 GeV, while the opposite-side lepton satisfies pt > 2.0 GeV. In 
addition, vertexing algorithms with efficiencies that decrease near the primary 
vertex increase the fraction of the longer-lived sample components. 

Figure |^ shows the like-sign fraction vs. proper time for the ALEPH analysis. 
Table |8| lists the results of fits to the oscillation frequency Am^, along with the 
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source and fractional variation of the parameters which contribute most to the 
systematic error. The and cascade fractions are important since they dominate 
the mistag rate. The ratio of to B^ lifetime is also important, since this can 
enhance or suppress the fraction of B~^ to B^ at long hfetimes, altering the shape 
of the hke-sign fraction and distorting the Am measurement. 



Table 8: 


Fit results and dominant systematic error 


for dilepton analyses. 


Experiment 


Dominant Systematic 


6{Amd) (ps-^) 


Anid (ps~^) 


ALEPH 


cascade frac. (15%) 
Bg frac. [y/o) 


±0.025 
±0.019 


0.452 ± 0.039 ± 0.044 




cascade frac. (15%) 
Tb, sample comp. 


±0.048 
±0.043 


0.503 ± 0.064 ± 0.071 


CDF(e/i) 


cascade frac. (15%) 
6-Baryon tract. (40%) 


±0.030 
±0.021 


0.450 ± 0.045 ± 0.051 


CDF(SLT) 


cascade frac. (25%) 
tb+/tbo{5%) 
charm dilut. (100%) 


±0.004 
±0.021 
±0.032 


0.500 ± 0.052 ± 0.043 
(includes Jet Q tags) 


DELPHI 


cascade frac. (6%) 
Tb+/tbo(4%) 


±0.032 
±0.017 


0.480 ± 0.040 ± 0.051 


L3 


cascade frac. (15%) 
^0 frac. (14%) 


±0.022 
±0.015 


0.458 ± 0.046 ± 0.032 


L3(IP) 


cascade frac. (15%) 
S° frac. (14%) 


±0.037 
±0.026 


0.472 ± 0.049 ± 0.053 


OPAL 


cascade frac. (15%) 
rB+/TBo(6%) 


±0.011 
±0.023 


0.430 ± 0.043 ± 0.030 



4-3 Lepton- Jet- Charge Analyses 

This section introduces the second commonly used opposite-side production flavor 
tag, the jet charge, as it is used in events with a lepton-tagged decay vertex. 

4.3.1 EVENT SELECTION AND PROPER TIME MEASUREMENT 

Events with a single lepton passing each experiment's standard lepton selection 
are considered. The lepton identification procedure follows that of the dilepton 
analyses for ALEPH DELPHI @, L3 |l| and OPAL [||, except that 



ALEPH raises the momentum cut on the electron to 3 GeV to match that of the 



muon. The CDF [78| data sample is the same as that of the lepton-SLT analysis 
described in the previous section. These experiments search for a secondary 
vertex in the lepton hemisphere and convert to proper time, as described above 
for the dilepton analyses. 
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SLD select leptons with p > 2.0 GeV and pf^ > 0.4 GeV (pp' > 0.8 GeV) for 
the inclusive lepton |84] and lepton-Z) [^] samples, defined below. 

SLD's vertexing algorithm involves combining all the (well-measured) tracks in 
a hemisphere. The position of the b vertex relative to the primary is estimated 

by L = where Xi is the vector from the primary vertex to the point on 

the lepton track that is closest to the track of particle i, and Wi is a weighting 
function that is the product of three weights. The first function suppresses tracks 
which are consistent with the primary vertex; it depends upon the impact pa- 
rameter (three-dimensional) of the track. The function approaches zero for small 
impact parameters and becomes roughly constant for large impact parameters. 
The second function suppresses tracks that have a poorly measured intersection 
with the lepton, and the third gives more significance to tracks for which this 
intersection is close to the intersection of the lepton track and jet direction. This 
method results in a vertex resolution with a core distribution of width 170^m. 

SLD uses a second method in which a D-like vertex is searched for inclusively 
in the lepton hemisphere. Tracks are first classified as non-primary if they have 
an impact parameter (three-dimensional) to the primary vertex of > 3.5a and a 
momentum greater than 0.8 GeV. The non-primary tracks (excluding the lepton, 
which is presumed to come from the B decay) are vertexed, and the resultant 
pseudotrack defined by the D vertex and momentum is intersected with the 
lepton. If this system forms a valid vertex, an attempt is made to add a primary 
track first to the B then D vertex. If this whole procedure fails, SLD searches 
for a single non-primary track that forms a valid vertex with the lepton. If this 
succeeds, the remaining tracks are vertexed to form the D as above. Although this 
is not an exclusive D reconstruction, we denote this analysis by ID to distinguish 
it from the above method. 

Apart from the typical vertex quality requirements, SLD also requires that the 
total charge of the tracks in both the D and B vertices be -1, or 1. Because 
of the high probability of correctly assigning tracks to secondary and tertiary 
vertices from the excellent tracking detectors, the SLD group can use the total 
charge at the B vertex to select B^ events with a high efficiency, while excluding 
a large fraction of the B~^ events. They estimate the sample contains 60.8% B^ 
and only 19.6% B+. 

SLD obtains the B momentum for the / hemisphere as follows. An attempt is 
made to include charged tracks only from the B decay by ordering the tracks in 
the jet based on the significance of their impact parameter to the primary vertex 
(in three dimensions). These tracks are combined until the invariant mass of 
the combination exceeds 2.0 GeV. This set of tracks, along with the lepton, give 
the charged particle contribution to the B momentum. The contribution from 
the neutrino is given by the difference between the beam energy and the total 
visible energy in the lepton hemisphere. Neutral energy is added based on the 
parameterization E^o = 0.7 Eq + O.OIE'q, where E^o is the neutral energy to be 
assigned to the B meson, and Eq is the total visible neutral energy in the b jet. 
This method results in a core resolution of 8% on the B momentum. 

The SLD /D-jet charge analysis fits to the decay length distribution directly, 
and so does not use event-by-event momentum information. 

The event selection is summarized in Table 0. 
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Table 9: Jet charge tag definitions for the inclusive lepton analyses, along with 
the right-sign probability parameterization as a function of Qtag • The mean value 
of Pr is quoted for rough comparison. 



Experiment Qtag Tag Defn. 



Pr{Q tug) 

(or Event Weight) 



Events <Pr> 



ALEPH fl' 

CDF @ 
DELPHI 

L3 g 

OPAL II 
SLD H] 



SLD 11] (ID) 



w 



1+Nd\Q 



IQtagI 
2 

I -^max 



62,320 0.53 



Qy^(O) - 10Q°PP(1) 



■>opp 



<Pr>, (|Qtag|>0.1) 
<Pr>, (|gtag|>0.12) 

From MC 
1 



1 + g-0.26Qtag 



1 + e 



-0.32Qtag 



120,700 
60,381 
8,707 
94,843 

2,609 
584 



0.56 
0.69 
0.72 

0.68 
0.68 



4.3.2 JET-CHARGE FLAVOR TAG 

The basic building block for the jet charge tag Qtag is either the weighted hemi- 
sphere or 6-jet charge, commonly defined as 

Qh,j{k) = — Pi,i = Pi-a 

Here a is the unit vector that is used to divide the event into two hemispheres 
(thrust, sphericity or &-jet axis), pi^i is the i*^ track's momentum component 
along this axis, and the sum runs over all charged tracks in the hemisphere 
{Qh, used by ALEPH, DELPHI, L3 and SLD) or jet {Qj used by CDF, OPAL). 
OPAL normalizes to the beam energy El^^^^ rather than J2iPi,i'i SLD has no 
normalization in the denominator; and L3 uses C{(f))qi in place of qi, where C{(p) 
is a function that deweights tracks that pass near an anode wire of the tracking 
chamber, for which the charge is not well determined, n is chosen to optimize 
the power of Qh to distinguish b from b quarks. For all experiments except SLD, 
— 1 < Qh < 1- Qni^^ = 0) is simply the mean track charge in the hemisphere, 
J27=i1i/''^ (sum of charges for SLD), and so is typically small. Qh{k = oo) is 
simply the charge of the track with the largest longitudinal momentum fraction, 
so <5j/(oo) = ±1. 

Qh is the basic building block for the jet-charge Qtag; however several experi- 
ments use a linear combination of the two hemisphere (jet) charges (with different 
K weights) as the tag variable. While jet charge is basically an opposite-side tag, 
the fragmentation tracks on the (possibly mixed) i?-vertex side still contain in- 
formation about the 6-quark flavor at its production, although the correlation is 
weaker than on the opposite side. Since one is searching for mixing on this side, 
it is desirable to find a combination that only depends weakly on whether the 
B meson has mixed (otherwise the right-tag probability would depend upon the 
mixed/unmixed state of the B, complicating the fit). The sum of the charges 
in the vertex hemisphere, (i.e. Q^^{k = 0)) is relatively insensitive to whether 
the B meson mixed, as only neutral mesons can mix. L3, for example, finds that 
the combination Qtag = Q°h^{K' = 0.4) — Q^^{k = 0) has maximal analyzing 



B MIXING 



37 



power. The combinations used by other experiments in the l-Qtag analyses are 
given in Table ^. SLD combines the jet-charge tag with their polarization tag 
described in Section |3^ . Finally, note that for OPAL, Qtag has been defined with 
the opposite sign and a scale factor of approximately ten. Figure shows the 
measured jet-charge distribution from SLD, along with the underlying b and b 
predictions from Monte Carlo studies. 




Pe signed cosO Qopp 

Polarisation Tag Jet Charge Tog 

Figure 10: 6-quark production flavor tags for the SLD experiment. Left: Polar- 
ization tag, based on the polar angle of the event axis. Right: Jet charge. The 
dashed lines are the distributions for b and b quarks, and the points are the data. 



One of the advantages of the jet-charge tag is its high efficiency. At LEP 
and SLD, the b jets are mostly produced back-to-back. Since the vertex-side 
jet must have tracks passing through the silicon detectors, it is well within the 
solid-angle covered by each experiment, and the opposite-side b jet is similarly 
well contained. A comparison between the /-/ and l-Qtag samples shows the latter 
contain approximately ten times more events. 

At CDF, however, the initial bb pair is not produced at rest, but with a po- 
tentially large momentum along the beam axis. A typical rapidity separation 
between the two b jets is Ar] ~ 1. The tracking chambers cover the region 
— 1 < 77 < 1, so quite often the vertexed b jet is contained within the acceptance 
of the tracking system while the other 6 jet is not. Even worse, another, non-6 
jet may be present, diluting the tag's effectiveness. These effects result in the 
average right-tag probability of 56%, rather than ~ 67% that is typical for the 
LEP and SLD experiments. However, a search for a secondary vertex in the non- 
lepton jet is performed. The subset of events with such a vertex has a greatly 
enhanced right-tag probability, namely 63%. 

4.3.3 FIT METHOD AND RESULTS 

The sign of the tag variable Qtag is used to categorize the sign of the 6-quark 
charge in the opposite hemisphere. The analyses are based on the time depen- 
dence of qc{t) = sign(Qtag) • sign{qi{t)). CDF, DELPHI, L3, OPAL and SLD 
perform likelihood fits to the same- and opposite-sign data, while ALEPH per- 
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forms a binned-x^ fit. The quantity < —Qh • Qiit) > is formed in 0.2 ps bins of 
proper time. This is simply the standard charge correlation function qdt) with 
each event weighted by | Qh \ ■ Although weighting by the dilution T> = 2Pr — 1 is 
optimal, note that CDF finds that for their procedure, in fact, T) ~ Qh- 

Events with |Qtag| near zero clearly have little power to separate opposite 
hemisphere b from b quarks. DELPHI and L3 simply cut out events with IQtagI < 
0.1 and IQtagI < 0.12 respectively. 

OPAL and SLD, on the other hand, parameterize PriQtag) based on Monte 
Carlo studies. For example, the separation between Qtag for b and b jets at 
SLD is shown in Figure [l^, which can be combined to give the b probability 
as a function of measured Qtag. CDF prefers to measure the shape of Pr{Qtag) 
from data. They measure the raw dilution, Praw = {^ss — J^os)/{^ss + ^^os), 
where SS refers to same-sign and OS refers to opposite-sign lepton and jet-charge. 
They find that the raw dilution can be parameterized by a simple straight line, 
-Draw(IQtagl) = «|Qtag|. Since the raw dilution is lower than the true dilution 
due to mixing and cascade semileptonic decays, the final fit includes an overall 
normalization, T> = N£)Drawi\Qtag\)- The normalization is N^i ~ 1.5 — 2.0, 
depending on the subsample. 

In addition to estimating the production tag right-sign probability for each 
event, OPAL also estimates the sample composition fractions event-by-event 
based on the value of the output of the neural net that is used to select lep- 
ton candidates, as described in Section 

The quality of the vertex-side lepton tag depends on the parameters discussed 
in the dilepton section. The sample composition on the lepton side is very sim- 
ilar to the dilepton case. The right-sign probability for this decay-point tag is 
typically estimated from the Monte Carlo predictions for the and cascade 
fractions. The tag purity of the jet charge is extracted from the data (see Table ^ 
for a summary and the number of events in each analysis). The systematic errors 
are dominated by the same effects which distort the relative amount of mixed 
and unmixed events as a function of the proper decay time, namely the sample 
composition parameters and the B^-B^ lifetime ratio. The results of the fits are 
summarized in Figure 12 and Table 10. 



^.^ D* — I and D* — Jet Charge Analyses 

As described in Section |3]^, a reconstructed D*^ meson can be used as a very 
pure tag of the charge of the b quark at its decay point; moreover, events with 
D*~^ are also enriched in B^ content. Further, the tracks forming the D*~^ can 
be vertexed and used to estimate the B decay length. 

The experiments search for the following decay modes: D*~^ — > D^t:^, with 
K^TT^ , K^ir'^Tr^ , or K^tt^tt^tt^ . All experiments exploit the small Q- 
value of the D*~^ decay: Md*+ — Mj^o — = 5.8 MeV, meaning the decay prod- 
ucts are basically at rest in the D*^ rest frame, which greatly suppresses the com- 
binatorial background. The quantity AM = M(all found decay products including 
M(all found decay products excluding vr^) is very insensitive to the momentum 
resolution of an individual track, since resolution effects mostly cancel in the dif- 
ference. This is true even to the extent that the vr*^ can be missing and a clean 
D* signal extracted. 

ALEPH JT^ fully reconstructs the three modes listed, including the vr" — > 
77. CDF [86[| reconstructs the Z)*+ only in the — > K~tt~^ mode, but adds in 
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Table 10: Fit results and dominant systematic errors for the lepton-jet charge 
analyses. 



Experiment 


Dominant Systematic c 


)(A?7i^) (ps ^) 


Am/7 fos ^) 


ALEPH 


rBo(3%) 


±0.012 
±0.016 


0.404 ± 0.045 ± 0.027 


CDF 


JD ' / ID \ ^ / 

right-sign frac, fakes (2%) 


0.021 
±0.014 


0.500 ± 0.052 ± 0.043 
(includes / tags) 


DELPHI 


fake frac. (20%) 
5-Baryon frac. (33%) 


±0.010 
±0.010 


0.493 ± 0.042 ± 0.027 


L3 


<Pr> (2%) 


±0.037 


0.437 ± 0.043 ± 0.044 


OPAL 


5° fraction (17%) 
rA,/rBo(8%) 


±0.011 
±0.007 


0.444 ± 0.029 ± 0.020 


SLD(0 


5° frac. (26%) 


±0.014 


0.520 ± 0.072 ± 0.035 


5-Baryon frac. (42%) 


±0.015 


SLD(/L') 


Fit binning and range 
Tracking efficiency 


±0.032 
±0.017 


0.452 ± 0.074 ± 0.049 



the — > K^Tx^Ti^ decay. DELPHI does not reconstruct the tt*^, resulting 
in larger background. OPAL [S7| finds the signal/background for the i^7r(7r'^) 
and K?>Tx modes to be too high, and excludes them from the analysis. Table ^ 
lists the modes reconstructed, as well as the signal-to-background ratio for each 
analysis. 



4.4.1 FLAVOR TAGGING 

The vertex-side h quark flavor at its decay point has the opposite sign of the 
reconstructed D*^ . The production flavor tags are either the opposite-side lepton, 
or opposite-side jet-charge tags which were introduced above. 



4.4.2 DECAY LENGTH AND PROPER TIME DETERMINATION 

These analyses have the advantage of a built-in decay vertex, namely that given 
by vertexing the D'^ decay products. Unfortunately, intersecting this with the 
7r+ from the D*^ decay to obtain the B decay point is not viable; the pion is 
emitted almost at rest in the D*^ frame, and so its direction is the same as the 
direction in the laboratory frame. All analyses begin by vertexing the D 
decay tracks. 

Only OPAL attempts to explicitly find the B decay vertex, as follows. After 
vertexing the tracks forming the D*"*" and creating a pseudotrack passing through 
this vertex with the D*'^ momentum, OPAL attempts to combine this pseudo- 
track with other tracks in the hemisphere that are likely to have come from the 
B decay. These tracks are selected based on their momentum and angle relative 
to the D*^ pseudotrack direction, and they must have an intersection point with 
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Table 11: Overview of the decay modes, sample sizes, flavor tag and mean right- 
sign tag probability for the D*-l and D*-jet charge analyses. 



Experiment D / Modes Sig/Bgnd 



Prod. Tag 



1 






840/565 






7r+ 






1555/1096 


+ 4- — 




7r+ 


358/520 




460/510 


7r+ 


1554/745 


7r+(^0) 


1370/1121 




1288/3904 


7r+ 


253/95 



<Pr> 



ALEPH 17|] 
ALEPH 17| 



CDF 



OPAL 13] K-TT+ 



I 

{p > 3.0 GeV 
> 0.75 GeV) 

Jet Q 
Qh{0.5)- 0.08 j:qi 



tag 



> 0.1 



I 

Pt > 8.0 GeV 
pf^ > 1.5 GeV 

Jet 

I 

Neural net* 



0.79 



0.76 



0.79 



0.65t 



0.79 



^Includes wrong-sign D* mistags. "''See Table * Approximately equivalent to 
cuts of p > 3.0 GeV and pf^ > 0.75 GeV. 



the pseudotrack that is consistent with a b decay. The B momentum is estimated 
using the same technique as in the OPAL dilepton analysis. 

CDF estimates the most likely B decay position based on its distance to the 
D vertex, the momentum of the D, and knowledge of the D and B lifetimes. 
The B momentum is estimated from the momentum of the D meson, using a 
correspondence determined from Monte Carlo studies. 

ALEPH and DELPHI fit the charge correlation to variables described below. 

4.4.3 FIT PROCEDURE AND RESULTS 

ALEPH fits directly to the charge correlation of the D* sign and opposite hemi- 
sphere sign (/ or jet charge) as a function of the decay distance, qdLDo). 
This involves combining both B and D decay time and momentum distributions 
in Monte Carlo studies to generate the predicted decay length distributions used 
in the likelihood fit. 

DELPHI fits to the charge correlation as a function of the sum of the proper 
decay times of the B and D meson. For /5=flight distance of the B, Z£)=flight 
distance of the D, and / ~ Ib + Id = total decay length to the D vertex, one 
defines a proper time t = ts + tn = ^^^B + ^J§'^D ~ This approximation 

is good to about a percent. DELPHI takes the average value Pb = 0.7£'beam to 
convert I to t. 

CDF fits to the mean charge correlation in eight bins of proper time while 
OPAL performs an unbinned likelihood fit as described previously. 

The vertex-side tag (the sign of the D^*^), has a very high purity - the main 
source of mistags comes from B — > DDg events (see Section ^) . The B^ fraction 
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in these samples is high, with a typical value of 83%. None of the experiments 
measure the production tag's right-sign probability event-by-event; instead they 
use the mean value < Pr >, which is left as a free parameter in the fits. The 



values of <Pr> are summarized in Table 11 



Figure 12 and Table 12 summarize the results of the fits. 



Table 12: Fit results and dominant systematic errors for the D*-l and D*-jet 
charge analyses. 



Experiment 


Dominant Systematic 


6{Amd) (ps-^) 


Anid (ps ^) 


ALEPH 
(combined) 


cascade frac. (15%) 
frac. (9%) 


±0.025 
±0.019 


0.482 ± 0.044 ± 0.024 


CDF 


B+ frac. (35%) 
Tc(25%) 


±0.037 
±0.021 


0.562 ± 0.068 ± 0.050 


DELPHI 


cascade frac. (6%) 
tb+/tbo(4%) 


±0.032 
±0.017 


0.523 ± 0.072 ± 0.043 


OPAL 


cascade frac. (15%) 
Tb+/tbo(6%) 


±0.011 
±0.023 


0.567 ±0.089 ±0.029 



4-5 D*l — I and D*l—Jet Charge Analyses 



CDF |8§, DELPHI [|g] and OPAL ^ perform analyses in which the vertex side 
of the event has both a D*^ and a lepton. A semileptonic B decay involving a 
D* has the lepton sign opposite the D* sign, e.g. D*^l~ , whereas combinatorial 
background contributes to both opposite-sign and same-sign. The production 
flavor tags used are the opposite-side lepton (CDF), or the opposite-side jet charge 
(DELPHI, OPAL). 

4.5.1 EVENT SELECTION 

CDF searches for D*^ candidates in the low-pt dilepton samples described pre- 
viously, using the same decay modes described in the above section and similar 
techniques. OPAL adds the — > K~ (tt^) decay. Neither experiment recon- 
structs the 7r°. 

DELPHI does not actually reconstruct the D*'^ at all; they only try to find the 
TT^ from the decay. Tracks likely to belong to the B decay are selected by the in- 
clusive vertex method in the lepton jet, exactly as described in the dilepton analy- 
sis in Section f4.2| . Just as in the case of the missing vr'^, forming the mass difference 
AM = M(A11 found decay products including vr^)— Af(All found decay products excluding 
allows a relatively pure D*^ signal to be extracted even if several B decay prod- 
ucts are missing. The larger background can be controlled via the extra handle 
of the charge correlation of the lepton with the tt^, . 



Table 13 summarizes the event selection. 
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4.5.2 DECAY LENGTH AND MOMENTUM MEASUREMENTS 

In D*l analyses, the lepton track, which is presumed to come from the B meson, 
can be used to reconstruct the B decay point. The decay products are ver- 
texed to reconstruct the D decay point, and a pseudotrack is formed with the 
D momentum passing through this point. This pseudotrack is intersected with 
the lepton trajectory to give the B decay position. OPAL solves for the decay 
length L using the position of the primary vertex, this secondary B vertex and 
a constraint based on the direction of the ID*~^ system. CDF performs a similar 
vertexing procedure in the r(/) plane and projects the decay length vector formed 
by the B and primary vertices onto the / — D* momentum direction to give L^-y. 

OPAL and CDF estimate the B boost as (3-/ = p^^*^ /Mb x /C(y^*^, M;£,*+), 
and Pjsind = ^ /Mb^ </C> respectively, where the correction factor, /C, is 
determined from Monte Carlo studies. The full distribution of /C is used to derive 
the momentum smearing needed for the likelihood functions. Typical resolutions 
on the decay time are 15-20%. 

DELPHI converts the decay length found by the inclusive vertex algorithm to 
proper time exactly as in the dilepton analysis. 



Table 13: Overview of the decay modes, sample sizes, flavor tag and mean right- 
sign tag probability for the D*l-l and D*l-]ei charge analyses. 



Experiment 




Modes 


Sig/Bgnd 


Production Tag 


<Pr> 


CDF [107| 


K- 
K- 
K- 


7r+ 

7r+(7r0) 


167/49 

190/226 

173/83 




0.67 


DELPHI II 


< 


-/ 


4132/1823 


Jet 


0.69 


OPAL 183] 


K- 


7r+ 

7r+(7r0) 


406/49 
794/225 


Jet 

(IQtagI > 1) 


0.72 



* See e/x and /i/i in Table 0. ^ See Table |. 



4.5.3 FIT METHOD AND RESULTS 

Each experiment performs a likelihood fit to the charge correlation of the charge 
of the vertex-side lepton and the opposite side lepton or jet-charge sign, as a 
function of the proper decay time. 

The decay-point tag is of course very pure, having both a D* and lepton require- 
ment. None of the experiments uses an event-by-event measure of the production 
tag's right-sign probability; instead, they use the mean value <Pr>, which is a 
free parameter in the fits. To increase this mean probability, OPAL and DELPHI 
remove low IQtagI events. The fitted values of <Pj.> are given in Table [l3| . 
Figure shows the DP— (mixed) fraction for the DELPHI 7^^,1 — 1 analysis. 



Figure O and Table 14 summarize the results of the fits. The systematic errors 



are dominated by the uncertainty in the B contamination. 
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Figure 11: Left: Fraction of vr,,,/ — Z events with like-sign (mixed) leptons, from 
the DELPHI analysis, and the fit results. Right: Right/wrong sign asymmetry 
from the CDF same-side tagging analysis and fit results. 



4-6 Other Inclusive Analyses 

There have been several attempts to increase the sample sizes beyond those in- 
volving leptons or D tags by using very inclusive methods, both in locating 
candidate B vertices and in flavor tagging. 

Both ALEPH [89 1 and SLD [^] have selected events using similar inclusive 
topological vertexing methods. SLD applies kaon and dipole final-state tags, 
using jet-charge and polarization for the initial tag, while ALEPH uses a highly 
efficient double jet-charge method. 



4.6.1 EVENT SELECTION 

The vertex algorithms used by ALEPH and SLD, rather than being driven by 
track intersections, search for viable secondary vertex positions. The ALEPH 
vertex algorithm is as described in the dilepton section, and is based on the 
difference of assigning all tracks to the primary vertex vs. assigning some to a 
secondary vertex. The SLD algorithm is similar; it searches for points in space 
from which several tracks originate with high probability. In both algorithms, no 
attempt is made to form separate B and D vertices. 

SLD finds a suitable vertex in ~ 50% of b hemispheres, compared to ~ 15% 
for charm and ~ 3% for light quarks. A vertex axis is formed by the direction 
between the primary and this secondary vertex. A final pass is made to assign 
tracks to the vertex based on their impact parameter to the vertex axis, and 
the distance along the vertex axis of the closest approach position. Finally, the 

quantity Mp^ = \J -|- pi + \pt\ is formed, where M is the mass of the tracks 
in the vertex ( assuming they are all pions) and pt is their total momentum 
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Table 14: Fit results and dominant systematic errors for the D*l analyses. 



Experiment Dominant Systematic 6{Am^) (ps ^) Am^ (ps ^) 



CDF 




frac. (30%) 

^ L>* frac. (50%) 


±0.020 
±0.027 


0.516 ±0.099 ±0.035 


DELPHI 




frac. (25%) 


±0.012 


0.499 ±0.053 ±0.015 


OPAL 




frac. (50%) 


±0.019 


0.539 ± 0.060 ± 0.024 



transverse to the vertex axis. Selected events are required to have Mp^ < 2.0 
GeV, which rejects most of the remaining charm and uds events, leaving 16,803 
vertexes with an estimated b purity of 93%. 

ALEPH uses the vertex as a seed and accepts tracks with Rsig = Sp/ {Sp + 
Ss) > 0.7, where Sp^s is the significance of the impact parameter of the track to 
the primary and secondary vertices respectively. Hemispheres with a vertex 
probability > 1% are accepted - if both hemispheres have a valid vertex, the one 
with the highest probability is selected as the vertex side of the event. 



Table 15: Overview of kaon, dipole and double jet-charge analyses. 



Experiment 


Events 


Decay Tag 


Production Tag 


ALEPH HI 


423,169 


Jet Q 


Jet Qt 


SLD @(Kaon) 


5,694 


Kaon 


Jet Q ± Polarization^^ 


SLD p0|] (Dipole) 


3,291 


Dipole 


Jet Q ± Polarization^^ 



t See Table |. 



4.6.2 INITIAL- AND FINAL-STATE FLAVOR TAGGING 

SLD determines the b flavor of the initial state by using both the jet charge and 
polarization methods, as described above. 

For the final-state tag, two methods were developed. First, the standard b 
s decay chain suggests that inclusive kaons should dominantly have the same 
sign charge as the b quark. This correspondence is diluted hy B ^ DD decays 
and by strange quark pairs produced during fragmentation. 

Kaons are required to be separated by more than 2.4cj from the pion hypothesis 
and by 3.2(T from a proton. The final-state tag is simply the sum of the charges 
of all kaon candidates, and so Qk > would tag a 5°, for example. The 12% 
of events with ^ Qk = are discarded. The probability of correctly tagging the 
charge of the b at its decay is 77% with this method. 

The second method involves constructing a charge dipole from the B and D 
meson decay products. SLD requires the total charge of tracks associated with the 
vertex to be zero, to enhance the B^ content of the sample as described previously. 
The direction between the primary and secondary vertices defines the vertex axis. 
A dipole is formed by summing over positive and negative tracks: 

5q = ^ , 
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where Li is the location of the distance of closest approach of track i to the 
vertex axis, and the weight Wi = siv? 6i/ aTi, where 6i is the angle of the track 
relative to the vertex axis and axi is the uncertainty on the impact parameter 
to the vertex axis. A —>■ l~^i'D~ decay, for example, will thus tend to have 
a negative dipole, whereas the equivalent B^ decay has a positive dipole. The 
probability of a correct tag is parameterized as a function of \Sq\, and reaches 
68% for large \5q\. 

ALEPH tags both hemispheres using the jet-charge algorithm previously de- 
scribed, using Qh{k = 0.5) for the opposite side hemisphere, and Qh{i^ = 1-0) 
for the vertex side. 



4.6.3 FIT METHOD AND RESULTS 

Both SLD analyses perform fits to the mixed/total (i.e. N^p^/N) fraction 
of the data, binned in decay length. Each event's contribution is weighted by its 
estimated dilution to maximize the sensitivity (see Section |3.3| ) . 

ALEPH converts the decay length to proper time using an event-by-event esti- 
mate of the momentum constructed from a weighted sum of the momentum of the 
charged tracks, energy from neutral clusters associated with the b jet (projected 
onto the vertex direction) and a correction due to the missing neutrino. They 
then fit to the weighted charge correlation < —Q°jf^{0.5)Q^^{l.O) > (t), i.e. the 
usual correlation function given by the signs of the hemisphere jet charges, with 
an event weight given by the product of their magnitudes similarly to the I- jet 
charge analysis described previously. 



The results are summarized in Table 16 and Figure 12 



Table 16: Fit results and dominant systematic error for double jet charge, dipole 



and kaon tag 


analyses. 






Experiment 


Dominant Systematic 


5{Amd) (ps"^) 


Arud (ps ^) 


ALEPH 


Frag, params. 
Afe frac. 


±0.015 
±0.012 


0.441 ± 0.026 ± 0.029 


SLD (Kaon) 


Fit method 
Kaon ID efficiency 


±0.042 
±0.035 


0.580 ± 0.066 ± 0.075 


SLD(Dipole) 


Fit method 
MC statistics 


±0.016 
±0.021 


0.561 ±0.078 ±0.039 



4-7 Same-side Tag Analyses 

As mentioned above, in an experiment at a hadron collider such as CDF, it is 
quite common that one b jet falls within the tracking acceptance while the other 
b jet does not. Hence a flavor tagging method which only relies on finding one b 
jet is attractive. Such methods are generally called same-side tags (SST). 
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4.7.1 EVENT SELECTION 

The only b samples available to CDF involve semileptonic b decays, selected by 
either the low-pt dilepton or high-pt single lepton triggers described above. Since 
this analysis is designed to only require one b, the single lepton samples are chosen 
as the basis. The SST algorithm used by CDF |91] requires a relatively pure 



meson sample, as discussed below. To obtain this, a reconstructed D*^ or Z)+ as- 
sociated with the lepton is required. Three decay modes, K'tt'^ , K^ir^(TT^), 
and K^'ir^TT^TT~ , (with the vr^ not reconstructed), are used, and the D*~^ re- 



construction follows the standard methods described in Section 4.4. The is 
reconstructed in the D'^ K~tt~^'k~^ mode. In addition to the D^*^+l~uX 
mode, the B~ — > D^l^vX mode, with K~^it~, is searched for, where the 

is required to not come from a D*~^. This signal is used to help estimate the 
amount of B~^ contamination in the (mostly) B^ sample. 



4.7.2 PRODUCTION FLAVOR TAG 

As described by the schematic picture in Figure ^, the SST tag used by CDF 
attempts to find the "first" charged pion in the fragmentation chain of the b 
quark. If the b quark forms a B'^, the first pion should be a vr"; conversely, a vr"*" 
would come from a b quark. However, if the b quark forms a B~ , the pion should 
be a TT~^, with a tt~ associated with a b, exactly the opposite correlation to the 
B^ case. Hence the charge of the pion, if it can be identified, determines the sign 
of the b quark, as long as it is known whether the decaying meson is a B^ / B^ or 

The SST pion is a fragmentation track, and so should originate at the primary 
vertex, not the B decay vertex. SST candidates are thus searched for near the 
B direction, defined as the / + D direction, with impact parameter significance 
to the primary vertex of less than three. String fragmentation models suggest 
that the fragmentation particles have little momentum transverse to the 6-quark 
direction, and so the track with the smallest pp^ is chosen, where pp' is the track's 
momentum orthogonal to the direction given by the I + D+SST pion system. A 
suitable candidate can be found in approximately 70% of the events. 



Table 17: Event selection and same-side tag properties from CDF. 







B^ sample, I'D 


) 


B sample, / 






K~7T+ 






Modes/Events 




^ i^-7r+(7r0) 


1754 
2515 


K'-K^ 2928 
(not from D*^) 








1997 




Tag Eflttc. 




70% 




70% 


Tag <Pr> 




0.59 




0.64 
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4.7.3 FIT METHOD AND RESULTS 

The decay length and conversion to proper decay time is performed as in the 



D*l — I analysis described in Section 4.5. 



The sign of the lepton in the ID^*^ decay tags the b quark charge at its decay 
point. For the predominantly sample, given by the ID*~^ events, the "un- 
mixed" DP+ sample has same-sign, /~ — Vggj^ and — '^ggj' events, whereas 
the DP— sample has opposite-sign events. For the -B"*" sample the assignment is 
reversed. 

Each sample is divided into six bins in proper decay length, and the mean value 
of the charge-correlation function < qdct) >= [N£,p^{ct)—N£,p_{ct)]/[N£,p^{ct)+ 
N£)p^{ct)] is evaluated for each bin. If the samples were pure B° and B+ , the 
first would show the usual cos Arudt dependence whereas the second would be 
flat. However, there is some cross contamination between the samples. First, B^ 
events with — > K~it^ decays end up in the B~^ sample if the tt^ from the 
D*^ decay is missed. CDF determined that this is the case for 15 it 7% of these 
B^ events. The most serious cross contamination, however, comes from decays 
involving the P-wave D** resonances (or non-resonant D^*^tt pairs), as shown in 
Figure 0. For example, the decay B~ D**^l^v followed by D**^ D*+tt^^ 
(denoting the vr from the D** as vr^,,,,) is classified as a B^ decay. A further 
complication is that the vr*,,, may be selected as the SST pion. As this example 
shows, this pion's charge is always correlated with that of the lepton such that 
the event is classified as DP+, which biases the charge-correlation function. The 
TT^^ particles, however, originate from the 5-decay vertex, whereas the SST pions 
are selected to come from the primary vertex, which suppresses this effect. 

The fraction of D** in semileptonic B decays, /**, and the branching fractions 
of the D** mesons are rather poorly known. CDF assumes f+f* + f** = 1, where 
/ and /* are the decay fractions into ID and ID*. CDF sets /** = 0.36 ± 0.12 
p2C and /*// = 2.5 ±0.6 Q. Further, there are four D** states of differing spin- 
parity, some of which decay to D*7r and some to Dtt. The relative abundances 
of the D** states in the B decay, then, have an effect on the amount of cross 
contamination among the samples by changing the fraction involving the D* 
state. This fraction is left as a free parameter in the fit and is found to be 
0.3 ±0.3. 



Table 18: Fit result and dominant systematic errors for the same-side tag analysis. 



Experiment 


Dominant Systematic 


d{Amd) (ps ^) Amd (ps ^) 


CDF |91 


1 


D** fraction(30%) 

D** D* fraction (100%) 


±0.031 0.471 ± 0.078 ± 0.034 



4 ■ 8 Summary of Mixing Results 

As this exhaustive summary of methods and results shows, a large number of 
different techniques are used to extract the oscillation frequency of B^ mixing. 
Within each single technique, each experiment typically has unique methods of 
measuring the decay length and momentum of the B meson. Further, in evalu- 
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ating systematics, the range of variation of variables, even those common to all 
analyses, is often not consistent. 

Therefore, combining all these results must be done with care. Since there are 
so many measurements, the final error can quickly become dominated by system- 
atic effects, and the handling of common systematics is especially important. A 
B Oscillation Working Group was formed [94| to combine the measurements in 
as consistent a manner as possible. 

Their results, shown in Figure 12, give a best fit of Arrid = 0.484 it 0.015 ps^^, 
including the time- integrated ARGUS and CLEO measurements. 



b2 mixing 



s 



That B'^ mesons undergo oscillations is clear. The time-integrated probability 
that a B^ mixes, Xd, is measured by the ARGUS and CLEO experiments, which 
operate at e^e~ colliders tuned to the T(4S') resonance. The T(4S') is not suffi- 
ciently heavy to decay to B'^B^ pairs, and so the only contribution to a mixing 
signal is from B^ mesons. At LEP, SLD and CDF, however, both B^ and Bg are 
produced, so time-integrated mixing measurements determine x = fdXd + fsXs, 
where fd ~ 40% is the B^ fraction in the b sample, and fs ~ 11% is the B^ 
fraction. Since x = 0.172 ± 0.010 and x = 0.118 ± 0.006 [|, a large value for Xs 
is required, indicating that B^ mesons do indeed mix. 

The problem is that B^ mesons oscillate at a high frequency, so Xs ~ 0.5 and 
Xs has no power to resolve Anig near this limit. Unfortunately, no experiment 
has yet succeeded in directly measuring the frequency of B'^ oscillations; only 
lower limits have been determined. 

For several reasons, including the notorious difficulty in combining exlusion 
regions from several experiments, an alternate method of fitting for Am^ was 
proposed by Moser and Roussarie [95|. This method, commonly called the ampli- 



tude method, provides a consistent fitting procedure for all experiments, making 
it much easier to combine the results. Essentially this method entails searching 
for a peak in the power spectrum of the data as a function of frequency. The 
likelihood function for a given data sample is constructed in the usual way, by 
estimating the amounts of various components of the sample. Each component's 
expected proper-time distribution is formed based on its true distribution, convo- 
luted with the decay-length and momentum-resolution functions, and combined 
in the appropriate fractions. The amplitude method introduces one difference in 
the construction of the likelihood function: the true probability densities of the 
Bg and B^ mesons are replaced by 

P{t) = -je-^'\l±AcosAmst), (34) 

i.e. an amplitude A is introduced in front of the oscillation term. 

The new likelihood function is then maximized as a function of A alone, i.e. 
all parameters, such as sample composition parameters, B hadron lifetimes, and 
even Am^, are fixed. This gives a measure of ^( Arris) and an error cr^(Ams) 
for that value of Arris . This fit is repeated for many different values of Arris to 
sketch out the shape of the ^(Arris) curve. If Arrig = Am*''"'^, the fit should 
return A = 1, within its errors. For Arris far from its true value, A should 
fluctuate around 0, within its errors. 

A given Am^ can thus be excluded at greater than 95% CL if .A+ 1.645cj^ < 1, 
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ALEPH D /I, 



ALEPH l/Qjet 
(91^14) 

ALEPH 1/1 
(91-94) 

ALEPH vtx/Qjet 
(91-95 pffel) 

CDF DVSST 
(92-95) 

CDFu/u 
(92-95) 

CDF l/l,Qjet 
(9*^5) 

CDF d'1/1 
(92-95) 

CDF e/u 
(92-95 pre!) 

CDF D'''/1 
(92-95 prel) 

DELPHI l/Qjet 
(91^14) 

DELPHI Tc'l/Qjet 
(91^14) 

DELPHI 1/1 
(91-94) 

DELPHI D'/Qiet 
(91^4) 

L31/1 
(94-95) 

L3 l/Qjet 

(94^5) 

L3 1/1(IP) 
(94-95) 

OPAL 1/1 
(91-94) 

OPAL l/Qjet 
(9V94} 

OPAL D'l/Qjet 
(90^14) 

OPAL d'/1 
(90-94) 

SLD K/QJet+pol 
(93-95 prel) 

SLD dipole/Qjet+pol 
(93-95 prel) 

SLD l+D/QJet+pol 
(93-95 prel) 

SLD l/QJet+pol 
(93-95 prel) 

Average of 25 above 

CLEO+ARGUS 
(Xj measurements) 



World average 




0.482 ±0.044 ±0.024 ps ' 
0.404 ±0.045 ±0.027 ps" ' 
0.452 ±0.039 ±0.044 ps ' 
0.441 ±0.026 ±0.029 ps ' 
0.471 tr„™ ±0.034 ps-' 
0.503 ±0.064 ±0.071 ps ' 
0.500 ±0.052 ±0.043 ps ' 

0.516 ±0.099 t^oisPs' 
0.450 ±0.045 ±0.051 ps ' 

0.562 ±0.068 t^otliPs"' 
0.493 ±0.042 ±0.027 ps ' 
0.499 ±0.053 ±0.015 ps ' 
0.480 ±0.040 ±0.051 ps ' 
0.523 ±0.072 ±0.043 ps ' 
0.458 ±0.046 ±0.032 ps ' 
0.437 ±0.043 ±0.044 ps ' 
0.472 ±0.049 ±0.053 ps ' 

0.430 ±0.043 tSps"' 
0.444 ±0.029 7„f,ps-' 
0.539 ±0.060 ±0.024 ps ' 
0.567 ±0.089 t"„'^j3ps"' 
0.580 ±0.066 ±0.075 ps ' 
0.561 ±0.078 ±0.039 ps ' 
0.452 ±0.074 ±0.049 ps ' 
0.520 ±0.072 ±0.035 ps"' 

0.483 ±0.016 ps ' 
0.483 ±0.033 ps ' 

0.484 ±0.015 ps ' 



B Oscillations 
Working Group 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 
Am . (ps"^) 



Figure 12: Summary of mixing results from the B Oscillation Working Group 



assuming a represents a gaussian error. An experiment's 95% CL lower limit on 
Anis is defined as the largest Am^ such that all smaller values have measurement 
probabilities below 5%. 

Of course, as with all limits, for a given error a^, fluctuations in the central 
value A can result in more or less favourable exlusion regions. The sensitivity of 
an experiment is therefore defined as the largest Am^ excluded if ^ = for all 
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Anis, i.e. the point at which 1.645cj^(Am<j) = 1. 

The lower hmits on Am^ and sensitivities of the analyses presented here are 
dominated by their statistical errors. The values quoted include a degradation 
due to systematic effects, but since these change the result very little, the major 
systematic errors are not listed here. 

Searches for B*^ mixing fall into two rough categories: high statistics but low 
-B^-purity analyses, based primarily on inclusive lepton samples, with typical B'^ 
purities of ~ 10%; and low statistics, high-purity (~ 60%) samples in which the 
Dg daughter of the B^ has been reconstructed. In order to obtain the maximal 
sensitivity to Arris, these experiments combine several production flavor tags. 



5.1 Inclusive Analyses 



These analyses closely follow the form of their B^ counterpart - the inclusive /-jet 



charge analyses described in Section 4.3 



5.1.1 EVENT SELECTION AND VERTEXING 
ALEPH ra, DELPHI [pTl, OPAL B 



and SLD [99, |100[ (with two analyses, 
referred to here as /-track and ID) search for Bg oscillations in the inclusive 
lepton samples described in Section [4.3| . A third SLD analysis (referred to as 
"dipole") is based on an inclusive vertex sample similar to the one used for the 
dipole-based B^ analysis described in Section 



Table 19 shows an overview of 



the event selection for these analyses. 



Table 19: Event selection for the inclusive-type B^ mixing measurements. All 
analyses except the SLD dipole are based on single lepton or dilepton samples. 



Experiment 



p{l) 
(GeV) 



WW 
(GeV) 



Events B^ fraction 




ALEPH |9e 
DELPHI 19/ 
OPAL i9l 
SLD 
SLD 



(/+track) 
(/+D) 
SLD pool (dipole) 



3.0 
3.0 
2.0 

1.0(e), 2.0(^) 
N/A 



1.25 
1.0 
0.7 
0.8 
0.9 
N/A 



33,023 
78,476 
53,050 
9,691 
2,009 
7,547 



10.4% 
10.5% 
10.5% 



15. 
15. 



ALEPH, DELPHI and SLD (/-track and ID analyses) search for a vertex by 
applying the same inclusive algorithms used in the B^ /-jet charge analyses, de- 
scribed in Section |]2| (ALEPH and DELPHI) and Section (SLD). For the 
SLD ID sample, a D-like vertex is required in the events. Demanding the total 
charge of tracks coming from the B and D vertices to be zero enhances the B'^ 
(and B^) content of this sample to 15.9%. DELPHI and SLD reconstruct the B 
decay length and momentum following the procedure used for the B^ analyses; 
in addition, SLD adds a few extra requirements designed to improve the decay 
length resolution and supress backgrounds from cascade 6 — > c — > / decays (such 
as the increased requirement). 

ALEPH adds requirements on the angle that the charm pseudotrack (produced 
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by the vertexing algorithm) makes with the lepton and with the jet, and on 
the mass of the B decay products. These cuts, which are not present in the 
B^ analysis, reduce the number of events by ~ 65%, but decrease the non-6 
background by a factor of ~ 4 and significantly increase both the decay-length 
and momentum resolution of the remaining events. The energy of the charm 
particle is estimated by clustering charged and neutral particles with energy > 0.5 
GeV (to reduce fragmentation clutter) with the charm vertex's tracks until the 
mass exceeds 2.7 GeV. The neutrino energy is estimated from the missing energy 
in the lepton hemisphere using the beam constraint as described above. The 

B momentum is then estimated hy ps = \J {Ec + + Ei)^ — m^. The core 
resolutions given by these methods are 280/im on the decay length and 7% on 
the fractional momentum resolution. 

OPAL uses a different vertexing algorithm than that used in the B^ inclu- 
sive lepton analysis. They attempt to find both the b and charm decay ver- 
tex simultaneously. They form a likelihood to which each track contributes 
{w/2) X (Pjj + Pc) + {1 — w) X Pp), where Pp, Pb and Pc are the probabilities that 
the track originates from the primary, b or charm vertex, as determined from the 
track's (three-dimensional) impact parameters to these vertex positions, and w is 
the probability that the track originates from the b vertex, determined from the 
track momentum and angle relative to the 6-jet direction. For the lepton track, 
w is set to 1. 

The inclusive vertexing algorithm used in SLD's dipole analysis differs from 
that used in the B^ analysis. In brief, it relies upon the fact that the D and B 
flight directions are almost identical. This implies that a single straight line (in 
three-dimensional space) can be found such that all the B and D decay tracks 
intersect it at one of two points. The algorithm first tries to find this line by 
minimizing a based on the intersections of all charged tracks with the line. 
Once this line is found, it is assigned a finite width and considered a "ghost" 
track. Tracks are again vertexed with this track or the primary vertex to build 
up the primary, secondary (B) and tertiary (D) vertices (see Reference |100| for 
the details of this method). The B momentum is estimated as in the ID analysis. 

5.1.2 FLAVOR TAGGING 

All the inclusive lepton analyses use the sign of the lepton to tag the b flavor 
at the decay point . The SLD dipole analysis defines a charge dipole 5Q = 
IbdX sign{QD — Qb), where Qb and Qd are the charges of the B and D vertices, 
and is the distance between them. As in the SLD B^ dipole analysis, 5Q > 
tags a B. 

ALEPH uses a combination of three production fiavor tags: 

• Opposite side lepton {p > 3.0 GeV) and opposite-side jet charge {QHiO-5)). 
See Sections and No cut is imposed on the lepton. 

• Fragmentation (SST) kaon: The kaon candidate, selected by dE/dx, is re- 
quired to be more consistent with originating from the primary vertex than 
the secondary vertex and have direction within 45° of the B^ direction. 

The data are categorized according to the presence of these tags. The lepton 
tag always takes precedence if it exists. If not, the kaon tag then used, if present; 
otherwise the sign of the jet-charge serves as the tag. The average right-sign 
probability of this tag is 71%. 
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DELPHI uses the following combination of production tag variables in the 
hemisphere opposite the lepton: 

• The opposite-side jet charge, Qh{0-^) (the hemisphere jet charge is defined 
in Section p!^ ); the opposite-side kaon charge, Qk (the jet charge evaluated 
for kaon candidates alone); and the opposite-side lepton, weighted by its pj*^' 
(see Section 

• The sum of charges of tracks not compatible with the primary event vertex 
and the sum of charges of tracks compatible with the primary event vertex 

In the lepton hemisphere, DELPHI uses: 

• The vertex-side jet charge: Qh{0-Q) 

• The fragmentation (SST) kaon or fragmentation (SST) (the rapidity of 
the highest-momentum kaon or candidate relative to the thrust axis that 
is compatible with being a fragmentation track, signed by its charge.) 

All this information, along with a tag based on the polar angle of the thrust 
axis (see Section p.2[ ), is combined into one variable. The average right-sign 
probability of this tag is 69%. 

OPAL also uses a combination of tags: 

• The opposite-side lepton - a lepton passing the standard OPAL require- 
ments. 

• The modified jet charge - a combination of jet charges from both the op- 
posite and vertex hemisphere, as well as fragmentation kaon information 
and opposite-side vertex charge. On the vertex side, Qvtx combines the 
unweighted 6-jet charge [Qh{0), with Qh given in Table ^ and the hemi- 
sphere charge [Qh{0) summed over all tracks in the hemisphere], as well as 
two other jet charges (with track weights given by neural nets which were 
fed the track's kaon and pion probabilities and momentum, direction and 
impact parameter information. 

The opposite hemisphere's tag, Qopp, combines the standard jet-charges 
Qni^-O) and Qh{0), the output of a neural net similar to that used in the 
vertex hemisphere, and, if available, the total charge of tracks in a viable 
secondary vertex. 

The tag variables Qvtx and Qopp are used to form the final variable 

Q 2(2vtxQopp -j^ 

QvtxQopp ~i~ (1 QvtxQopp) 

OPAL estimates this tag variable is 40% more eff'ective than the combined jet 
charge used in the /-jet charge analysis. 

All three SLD analyses use the same combination of jet-charge and polarization 
(polar angle) initial-state tags described in the /-jet charge section. 



5.1.3 FIT METHOD AND RESULTS 

The sample composition of the lepton hemisphere is determined event-by-event. 
It is similar to that of the inclusive lepton samples described in Section The 
resulting increase in sensitivity is equivalent to having approximately 30% more 
events. All experiments perform amplitude fits to the data. The individual lower 



limits on Am^ and analysis sensitivities are listed in Table 22 
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5.2 Analyses with Reconstructed Dg Decays 

The remaining analyses take a more exlusive approach, ranging from partially 
reconstructing the Dg decay in DJ decays, to fully reconstructing the 

Dg decay but partially reconstructing the Bg in Bg DglX, DghX decays, to 
fully reconstructing the Bg decay. Clearly the latter will give ideal vertex and 
momentum determinations, since all the Bg decay products are found; however 
the very small exclusive Bg branching ratios mean that the event sample is rather 
small. 

ALEPH and DELPHI reconstruct many Dg decay modes, summarized in Ta- 



ble 2C. The strange mesons are reconstructed in their usual accessible modes: 
(j) K+K~, K-TT+, K*^ K^Tr+, and the / resonance as /(980) ^ 

7r"'"7r~. For details on the selection criteria, see the references indicated in Ta- 
ble p. 



Table 20: Dg decay modes searched for by various experiments. 



ALEPH 


DELPHI 


CDF 


Dg mode 


Dgl 


Dgh 


Dgl/cpi 


Dgh 


Full B'^g 






Y 


Y 


Y 


Y 


Y 






Y 




Y 










Y 




Y 




Y 






Y 


Y 


Y 


Y 


Y 




K*^K*+ 


Y 




Y 




Y 




K^gK+ 


Y 


Y 


Y 




Y 




/(980)7r+ 
4>lui{l = e,iJ.) 


Y 


Y 


Y 




Y 






Y 


Y 



5.2.1 EVENT SELECTION 

ALEPH |0T|, |102| and DELPHI |To|, |lO|] reconstruct 5° Djl+X and 5° 



Dg h'^X, with the Dg decaying in the modes given in Table In addition, 
DELPHI searches for fully reconstructed Bg (up to a missing 7 or vr^) in the 
following modes: D^tt^ , Dj , D'^K^tt^ , and D^K^af, with af /3°7r+ 
7r"'"7r~7r"'" and D^ K^ir^ , K^tt^tt^tt^ . 

Note that for Bg — > D* Dg^, Dgir^ transitions, the photon or vr*^ may not 
be reconstructed. Thus the mass spectrum of the Bg candidates is expected to 
have a sharp peak at the Bg mass, and a wider satellite peak at lower mass from 
D'^ decays. DELPHI finds 8 it 4 signal events out of the 11 events in the Bg mass 
region, and 15 it 8 signal events out of 33 in the satellite region. 

In addition to the exclusive Dg modes listed, for events with a lepton, DELPHI 
also searches for a partially reconstructed Dg of the form Dg (ph'^X (which 
was not found as one of the above modes). 

CDF |105|] searches for cj) — > K~^K~ candidates in the low-pt dilepton samples 



(/u/x and e/x) described in Section 4.2. The lepton on the 4> side is required to have 
pp' > 1.0 GeV. A charged hadron near the pair is required, and the effective 
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Table 21: Summary of the Bg decay modes, sample sizes and purity, and right- 
sign tag probability for analyses with a partially or fully reconstructed Dg. 
Experiment Bg Mode Events Bg frac. <Pr> 



ALEPH |T0|, [T02| 



CDF [105| 



DELPHI 1T03|, [104 1 



Dgl 


277 


66% 


0.73 


Dgh 


1,620 


22% 


0.74 


(phi 


1,068 


61% 


0.76 


Dgl 


436 


53% 


0.78 


(phi 


441 


9% 


0.78 


Dgh 


2,953 


20% 


0.78 


B'g (full) 


44 


51% 


0.78 



masses are restricted to be 1.0 < mif,h < 2.0 and m^/jj < 5.0, consistent with 
decay kinematics for Bg Dglu and Dg — > (j)hX. 
Table || gives an overview of these data samples. 



5.2.2 DECAY LENGTH AND PROPER TIME 

A B decay vertex is easily obtained in these events. First the Dg decay products 
are vertexed and a pseudotrack formed, which passes through this vertex and 
has a momentum given by the Dg momentum. This track is intersected with the 
other track available, / or hadron, to give the B^ decay point. 

For the Dgl samples, ALEPH and DELPHI estimate the B momentum based 
on the momentum of the Dgl system and apply a correction for the neutrino ob- 
tained from the missing energy in the hemisphere. For the Dgh sample, ALEPH 
adds charged and neutral particles to the Dgh system and selects the most likely 
combination based on the effective mass. The momentum of the fully recon- 
structed decay is simply measured from all the decay products. CDF corrects 
the measured momentum of the (j)hl system for the missing energy using a mean 
value derived from Monte Carlo simulation. 



5.2.3 FLAVOR TAGGING 

All experiments use the sign of the lepton on the vertex side as the final-state 
flavor tag. ALEPH and DELPHI use the combination of production flavor tags 
described in the inclusive lepton section above, whereas CDF uses the sign of the 
opposite-side lepton. 



5.2.4 FIT METHOD AND RESULTS 

DELPHI estimates the right-tag probability event-by-event by using the shape 
of the Prix) distributions, which are determined from Monte Carlo simulation. 
This improves the effective mean right-tag probability, e.g., from 0.745 to 0.78 
for the 1994-1995 Dgl data. The effectiveness of the combined tags can be seen 
by comparing the right-sign tag probabilities in Table 21 to the individual tags 
used in the B^ analyses. 

All the experiments perform an amplitude style fit, with the Amg exclusion 



regions and sensitivities given in Table 22 
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Table 22: Excluded regions of Am^ and the sensitivity of each analysis. 



Experiment 



Excluded Amg 
region (95% CL) 



Sensitivity 



9.5 ps"^ 

6.7 ps-i 
4.1 ps-i 

5.1 ps-i 
6.5 ps-i 

8.2 ps-i 
3.2 ps-i 
7.2 ps-i 

3.8 ps-i 
3.5 ps-i 
5.4 ps-i 



ALEPH / [^1 



ALEPH Dsi 

ALEPH D.h^^m 
CDF (pi - / JlO| 
DELPHI / 1^ 
DELPHI Dsl |10| 
DELPHI _B0 + Dsh |0|] 
OPAL / 
SLD / 

SLD ID im 

SLD dipole pO 



< 1.3 ps-1,2.0 



< 9.5 ps"^ 

< 6.6 ps"i 
< 3.9 ps-i, 6.5 - i 

< 5.8 ps-i 

< 4.7 ps"i 

< 7.5 ps 

< 4.0 ps 

< 5.2 ps 
- 8.6 ps" 

< 5.2 ps 

< 5.2 ps 



ps" 



9.8 - 12.2 ps" 



Table ^ summarizes the excluded regions of Amg from each of the analyses 
presented, along with their sensitivities. The B Oscillation Working Group has 
combined these, and finds that Am^ > 14.6 ps"^ at the 95% CL |Q. Figure 
contains the combined amplitude plot showing the value of the amplitude A in 



Equation 34 versus assumed Amg, from which the lower limit and experimental 



sensitivity can be read. 
6 5? LIFETIME DIFFERENCE 



s 



Another, more indirect way of searching for B'^ mixing utilizes the relationship 
between the mass difference of the two eigenstates, Am^, and the decay rate 



difference AF, [|lO§ 



AT, ^ _3^^^ 
Am, 2 m?,,^255 

The QCD correction factors r] have recently been estimated, giving AFs/Am^ = 
(5.6 ±2.6) X 10-3 |0|. 



A large Am^ corresponds to a large lifetime difference between the heavy and 
light states. Even a moderately large Am^ ~ 20 would result in a 17% difference 
in their lifetimes. 

This difference has been searched for in two ways. First, B^ decays to final 
states that are not CP eigenstates should contain a mixture of both Bh and Bl 
components. Several experiments have fit their B^ lifetime distributions to a sum 
of two exponentials of decay widths Tg ± AFs/2. From their sample of ~ 600 
S° ^ l+DjuX events, CDF Q finds AfJf^ = 0.34t|]|^ (statistical only). 



corresponding to an upper limit of AFs/F^ < 0.83 (95% CL). DELPHI 
using their D^P and D^hP samples from the B^ mixing analysis (Section 
finds AF,/F, < 0.42 (95% CL). 



lOJ 
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^ 3 

s 

a 2.5 
E 
< 

2 



' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' 
World average (prel.) 

j data±la ▲ 95% CL limit 14.6 ps"' 
1.645 a > sensitivity 14.6 ps"' 

■ □ data ± 1.645 o 

I I data ± 1.645 o (stat only) 




8 10 12 14 16 18 20 



Am (ps' ) 



Figure 13: B Oscillation Working Group's |94| combined plot of fitted amplitude 
versus Am^, for the Bg mixing analyses presented in this section. 



L3 |10S| ] fits 445,000 events, selected using an inclusive vertexing algorithm, to 
a sum of B^, B'^, Ai, exponential decay distributions and a two-lifetime -B^ com- 
ponent, resulting in Ar^/Fs = O.OOI^q qq (statistical only). Including systematic 
effects, L3 estimates an upper limit of Ar^/Fs < 0.67 (95% CL). 

A second method to search for a lifetime difference is to measure the lifetime of 
Bg decays to modes that are CP eigenstates or are mixtures dominanted by one 
CP state. These decays arise primarily from either Bh or Bl decays, depending 
on the CP eigenvalue of the final state. The difference between the lifetime 
measured in this sample and the mean lifetime measured over a mixture of final 
states measures AF. For example, ALEPH ||110| ] found 32 it 17 decays of B^ 

_^ ^^x, from which they extracted t{Bs) = 1.42 ± 0.23 ± 0.16 ps. 
Assuming this is a pure CP-even decay, comparison to the mean B^ lifetime 
of TB, = 1.55 ps gives AF./F^ = 0.17 ± 0.35. 

CDF measures the lifetime using 58 it 12 decays of B*^ — > J/i^(j) finds 
t{Bs) = 1.34 lb 0.23 lb 0.05 ps [111]. This final state is not necessarily a pure CP 
eigenstate, and using a transversity analysis |112] CDF finds that the CP-odd 
component represents 23ib 19% of the sample. The higher statistics samples soon 
to be available at the Tevatron will clearly benefit this type of analysis. 

The DELPHI result, combined with the mean S° lifetime, gives an upper limit 
of Am, < 48 ps"^ (95% CL). 
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7 CONCLUSION 



Mixing between the particle and antiparticle states of the and B^ mesons has 
several consequences. It leads to physical states that are no longer eigenstates of 
flavor, particle-antiparticle oscillations can generate CP asymmetries between B 
and B decays and their frequency is relatively directly related to CKM matrix 
elements. The B system is unique in having two such particles readily available, 
which can be studied and compared. 

Several experiments have succeded in observing the time-dependence of B^ 
flavor oscillations, using a wide variety of techniques. The measurement of the 
B^ oscillation frequency will have to wait for the next generation of experiments 
- current experiments are only able to place a lower limit. The principle results 
presented here are: Am^ = 0.484 it 0.015 ps~^ and Am^ > 14.6 ps^^. 



With the results from Section [2. 
IbJIb. = 1.15 ± 0.06, BbJBb, 



/b, = 200 ± 16 MeV, Sg^o 
1.01 ± 0.01 and ^7(174)" 



0.55 lb 0.01, the CKM matrix element \Vtd\ and ratio |Vj 



rjB 

extracted according to Equation 24 



ts 



= 1.37 ±0.08, 
167 ± 5 GeV, 
\/\Vtd\ can be 



Arrid 



Anis 
Anid 



24.4 



200 MeV 



B 



NLO 



1.37 



±5.0] \Vtd\^ (eV) 



(1.36 ±0.14) 



V, 



ts 



Vtd 



corresponding to \Vtd\ = (3.6 ± 0.4) x lO^^ and \Vts/Vtd\ > 4.7 (95% CL). 

The future direction for B mixing measurements is clear. Measuring Anis, 
either directly through a time-dependent oscillation analysis, or through a heavy- 
light lifetime difference is crucial. The prospects for this occurring in the next few 
years are good, as the upgraded CDF and D0 detectors are expected to provide 
significantly improved statistics and sensitivity to measure Anis- 

It is equally important to make further progress in calculating fs and Bb- The 
lattice calculations of these constants are maturing rapidly, and with unquenched 
simulations becoming available, more precise estimates with a full evaluation of 
systematic uncertainties should be within reach. 

These measurements are an important part of exploring the unitarity of the 
CKM matrix, along with the various -B-sector CP violation searches which will 
take place over the next few years. The study of B mixing will surely remain 
vibrant for years to come. 
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